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BACKLASH  IN  SECON1    ORDER 
FEEDBACK  CONTROL  SYSTEMS 

1 .     Introduction 

Whenever  components  are  mechanically  connected  by  a  linkage 
or  a  gear  train,   there  exists  the  basic  problem  of  providing  smooth 
operation  without  looseness  in  the  connecting  parts  of  the  mechanical 
system.     In  general  this  basic  problem  is  not  completely  solved,   and 
the  connecting  parts  are  free  to  move  independently  over  a  limited 
range.     This  phenomenon  is  commonly  called  backlash   and  its  effects 
are  usually,,   but  not  always ,   undesirable.     In  particular ,  when  back- 
lash  occurs  in  a  feedback  control  system  it  often  causes  a  small 
amplitude  "dither"  or  limit  cycle,   and  when  the  system  reaches  static 
steady  state  without  a  limit  cycle  the  accuracy  of  the  system  is 
affected  by  the  backlash. 

In  the  usual  feedback  control  system,   the  backlash  occurs  in 
two  places:    in  the  power  drive  gear  train  between  motor  and  load, 
and  in  the  gear  train  between  load  and  output  measuring  device  .     It 
is  not  difficult  to  establish  the  pertinent  differential  equations  for 
either  of  these  cases,   nor  is  there  any  conceptual  difficulty  in  the  solu- 
tion of  these  equations,   but  considerable  labor  is  required  for  the 
solution  of  even  the  simplest  cases.     Furthermore,   the  solution  of 
one  or  two  specific  cases  does  not  provide  perspective  for  generaliza= 
tions ,   and  the  multiplicity  of  parameters  involved  makes  detailed 


analyses  unattractive.     This  report  deals  with  some  initial  results  in 
a  long  range  project  to  study  the  backlash  phenomena  quantitatively, 
with  the  ultimate  goal  of  providing  charts  which  will  be  useful  in  the 
engineering  analysis  and  design  of  feedback  control  systems  involving 
backlash  . 

In  general,   the  simplest  linear  feedback  control  system  has  a 
second  order  differential  equation,   and  is  commonly  referred  to  as  a 
second  order  system.     Backlash  may  and  does  occur  in  such  systems. 
While  it  is  obvious  that  most  feed)  ack  control  systems  are  third  order 
or  higher,   the  study  of  backlash  in  this  report  is  restricted  to  second 
order  systems  only,    so  as  to  provide  a  firm  foundation  for  future 
studies  while  limiting  the  variable  parameters  to  a  variety  which  makes 
numerical  calculations  feasible.     In  like  manner  the  mechanical  para- 
meters  are  restricted  to  inertia  (mass)  and  viscous  friction;  couloumb 
friction,   though  of  some  practical  importance,   is  not  considered  quanti- 
tatively and  is  permitted  in  a  qualitative  sense  in  only  one  special  case. 
It  is  recognized  that  the  mechanical  inertia  and  viscous  friction  may 
be  subdivided  and  appear  on  both  sides  of  the  backlash  in  varying  pro- 
portions,   such  cases  are  considered  in  detail.     It  is  also  recognized 
that  the  gear  tooth  surfaces  must  hit  each  other  after  the  meshing  gears 
have  gone  through  the  dead  space  (backlash),   and  this  impact  may  be 
either  plastic  or  elastic;  this  effect  is  also  studied  in  detail.     The 
mechanical  load  on  the  gear  train  is  considered  to  be  only  inertia,  and 


viscous  friction.     This  combination  contains  only  one  energy  storage 
device  (inertia)  and  is  therefore  a  first  order  load,     In  some  cases 
the  mechanical  connection,   e.g.    shaft,  may  act  as  a  spring,   thus  pro- 
viding  a  second  order  load  which  is  capable  of  resonance  if  excited  by 
periodic  motion  of  the  proper  frequency.     This  case  is  not  considered 
but  is  reserved  for  future  study. 

The  two  basic  tools  for  studying  the  nonlinear  effects  of  backlash 
are  the  phase  plane-phase  space  technique  and  the  describing  function 
method.     The  describing  function  method  has  been  applied  to  a  number 
of  backlash  studies,   and  it  is  intended  to  devote  some  future  efforts  to 
an  extension  of  such  studies;  however,   the  investigations  reported  on 
here  are  based  solely  on  the  phase  plane-phase  space  technique.     For 
second  order  systems  the  phase  space  is  two  dimensional,   that  is, 
it  reduces  to  a  phase  plane s   and  nonlinear  as  well  as  discontinuous 
systems  can  be  analyzed  using  the  welll       vn  graphical  methods  of 
integration  (method  of  isoclines,    Lienard' s  construction^   the  Delta 
method)  which  make  the  phase  plane  a  powerful  tool.     This  provides 
a  convenient  starting  point  for  the  analysis,   though  the  graphical 
methods  do  not  apply  to  higher  order  systems, 

The  phase  plane  technique  with  graphical  integration  was  used  in 
initial  studies  recorded  here  but  was  applied  to  only  the  simpl« 
cases.     The  results  obtained  and  the  procedures  developed  provided 
sufficient  insight  to  permit  formulation  of  a  rather  comprehensive 


digital  computer  program  to  study  the  effects  of  many  parameters  on 
the  response  of  a  given  system       The  majority  of  the  data  presented 
here  was  obtained  using  a  digital  computer  program;  however,   any  of 
the  results  so  obtained  could  be  duplicated  by  graphical  methods. 

The  digital  computer  was  programmed  to  provide  a  complete  solu- 
tion of  the  differential  equations,   and  sufficient  data  was  computed  to 
provide  a  complete  and  very  accurate  phase  trajectory.     However,   for 
most  of  this  initial  phase  the  question  of  primary  interest  was  the 
development  of  criteria  for  the  existence  of  limit  cycles       The  computer 
program  was  therefore  provided  with  a  subroutine  which  printed  out 
only  the  maximum  position  per  cycle.     This  was  done  to  decrease  the 
read  out  and  print  out  times.     By  a  simple  auxiliary  command,    complete 
data  for  the  phase  trajectory  was  made  available;  however,   the  com- 
plete phase  trajectory  was  solved  for  only  a  restricted  number  of  cases 
because  of  the  time  involved.     As  a  result,   the  limit  cycle  studies  are 
detailed  and  fairly  complete  (more  data  would  improve  the  results  but 
should  not  introduce  any  startling  features);  but  many  questions  con- 
cerning response  during  the  transient  period  have  not  been  completely 
answered.     Much  of  the  remaining  work  (within  the  restricted  scope 
of  this  initial  phase)  is  just  a  matter  of  data  reduction. 


2 .     THE  NO  LOAD  CASE 

For  some  systems  one  may  neglect  inertia  and  viscous  friction 
on  the  output  side  of  the  gear  train  containing  the  backlash.     This  is 
arbitrarily  called  the  No  Load  Case.     Two  practical  variations  are 
common.     First,   the  most  commonly  encountered,   is  the  case  of  the 
instrument  servo  which  is  driving  some  simple  indicating  or  measur- 
ing device  such  as  the  wiper  arm  on  a  potentiometer.     In  this  case 
the  actual  driven  load  has  a  negligibly  small  inertia  and  viscous  fric- 
tion.    There  is  usually  some  couloumb  friction  in  this  load,   but  it  is 
often  too  small  to  affect  the  differential  equations  though  large  enough 
to  prevent  jitter  in  the  low  inertia  driven  load.     The  second  case  of 
interest  here  is  thc-t  in  some  larger  systems  with  driven  load  of  appre- 
ciable inertia,   the  power  gear  train  has  negligible  backlash,   but  there 
is  a  measuring  device  geared  to  the  load,   and  there  is  backlash  in  the 
gear  train  between  load  and  measuring  device.     Note  that  while  these 
two  cases  arise  in  systems  which  are  physically  different,   they  are 
mathematically  equivalent  and  thus  one  analysis  suffices. 

The  block  diagram  description  of  the  system  is  given  in  Fig.    2-1. 
L  it  is  assumed  that  a  step  displacement  input  is  applied  to  the  system 
when  the  gears  are  in  contact  and  the   system  is  driving  linearly,   then 
the  basic  linear  differential  equations  are: 
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The  dynamic  response  of  the  system  is  shown  on  the  phase  portrait 

of  Fig.   2-2,   and  equations  2-1  and  2-2  describe  the  first  portion  of  the 

trajectory  from  MtoN.     At  point  N}   E  =  O  and  0    =  Os   the  motor  re- 

verses  and  runs  in  the  opposite  direction,   increasing  the  velocity  and 

displacement  of  the  motor  shaft.    However  the  output  shaft  cannot  move 

and  the  error,   E,   cannot  change  until  the  backlash  is  taken  up.    The 

operation  during  this  interval  is: 
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It  is  readily  seen  that  the  horizontal   or  E-axis  is  a  locus  of  points 
at  which  the  describing  differential  equation  changes  discontinuously , 
and  it  may  then  be  called  a  "dividing  line."     This  particular  dividing 
line  is  a  "separation  line"  i.e.   as  the  state  point  crosses  this  line  the 
physical  system  goes  into  the  backlash  region  and  the  drive  member 
is  mechanically  separated  from  the  load  member.     During  the  time 
interval  when  equations  2-4  and  2-5  apply,   and  the  backlash  is  being 
"taken- up";  the  load  member  does  not  move,   E  and  E  do  not  change. 
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but  0      and  0       do  change.     On  the  phase  plane  of  Fig.    2-2  the  vertical 

mm  r  ° 

line  from  N  to  P  indicates  the  first  interval  of  this  type       It  can  be 
argued  whether  this  line  is  or  is  not  validly  a  segment  of  the  phase  tra- 
jectory s  but  the  physical  interpretation  is  quite  clear:     E  and  E  are 
constant  at  point  N  until  the  backlash  is  completely  taken  up,   then  as 
the  gear  teeth  hit  (with  inelastic  impact)  the  load  instanteously  acquires 
a  velocity  and  the  state  point  jumps  dis continuously  from  N  to  P.     Once 
the  backlash  is  taken  up  the  system  once  more  operates  linearly  and 
the  segment  of  phase  trajectory  from  P  to  Q  is  again  governed  by 
equations  2-1  and  2-2.     The  point  P  is  one  point  on  another  dividing 
line,   this  is  the  line  at  which  the  backlash  is  taken  up  and  the  system 
physically  recombines,    i.e.  ,   the  mechanical  connection  is  restored. 
Thus  the  line  may  be  called  a  backlash  "recombination  line."     This 
line  can  be  constructed  accurately  only  by  a  point-by-point  calculation. 
However,    equation  2-5  can  be  rewritten: 

where  t  a    is  the  time  interval  required  to  take  up  the  backlash.     For 
t  a        00  the  slope  of  the  dividing  line  is  -1/2   y  .     Thus  for  various 
values  of  the  backlash,^  ,   the  recombination  line  may  be  represented 
on  the  phase  plane  as  shown  in  Fig.   2-3.     Only  the  tangent  line  is 
accurate  (unless  detailed  computations  are  used),  but  knowledge  of 
this  tangent  line  permits  derivation  of  a  stability  criterion. 
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FAMILY    OF    BACKLASH    DIVIDING    LINES   WITH 
COMMON    TANGENT   AT   THE  ORIGIN. 
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2 .  1    A  Stability  Criterion  for  the  No  Load  Case 

If  equations  2-1  and  2-2  are  combined  and  nondimensionalized 
forming  a  single  equation  with  error,   E,   as  the  independent  variable, 
this  equation  can  be  manipulated  and  integrated  to  define    he  linear 
phase  trajectory  on  the  phase  plane  as: 

E2  +  2^EE+u.n2E2  ^  _,      ^        E 

tan 
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Using  the  linear  transformation 
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Letting      Y=    ©     sin  -ty  ,    X=  p     cos  yty    ,    \|f     =  tan        Y/X 


p=De^1-S*  <2    10> 

The  backlash  dividing  lines  are  located  on  the  X  -   Y  plane  by  noting 
that 
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but  the  dividing  line  for  separation  is  the  E-axis,   for  which  E=Os   and 
for  this  line 
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The  backlash  dividing  line  for  recombination  is  not  readily  determined, 
but  the  tangent  to  this  curve    A  on  Fig.    2-3,   is  easily  transferred: 


max 
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y    2r2-  i  (2~14) 

and  for  \.  E  /  E  wn  =  —  1  / 2  (»  so      ~= =  tan  \h 

max  "  °  X  -  ^max 
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which  may  be  expressed  as  an  acute  angle 

tan    yh  =      1   ~~Z* (2-15) 
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2fVl_    tZ 

In  Fig.    2-4  the  backlash  line  for  recombination  is  sketched,   not  calcu- 
lated,  but  the  basic  condition  for  increasing  oscillations  is 
Pl    cos    y    l  =    p2    cos    >pt  (2l6) 

Assuming  that  the  angles  are  measured  from  the  negative  real  axis  in 
Fig     2-4  and  that  p       and  Z     are  known: 
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Substituting,  .  \U  yu 


-$(n-*T-Vl)  (2.20) 
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The  angle^      may  be  chosen  arbitrarily;   and  by  inspection  of  Fig.    2-4 
it  is  apparent  that  if  the  required  inequality  exists  it  must  be  most  pro- 
nounced when  the  point  (p    >^,)  is  located  near  the  origin  on  the  tangent 

line,   in  which  case  Vf   =  V  •     Substituting  this  value  for    Uf ,: 

'  1       'max  '  1 


(2-22) 
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Equation  2-22  has  been  plotted  in  Fig.   2-5,   and  the  stability  point  is 
marked.     This  shows  that  the  system  with  backlash,  but  no  load,   will  have 
a  limit  cycle  if  £  =    0.29,   and  the  existence  of  the  limit  cycle  is  inde- 
pendent of  the  amount  of  backlash,    so  long  as  some  backlash  is  present. 
Obviously  the  amplitude  of  the  limit  cycle  will  be  del  i.  rmined  by  the  amount 
of  backlash. 
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FIG.  2-5    GRAPH    OF    STABILITY    CRITERION    FOR    BACKLASH, 
k  s_e_-a  (II-  tan"'  b-  tan"'  a ) 
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3.     PHASE  PLANE  ANALYSIS  C  F  BACKLASH  IN  A  SECOND  ORDER 
SYSTEM  WITH  DIVIDED  LOAD  AND  PLASTIC  IMPACT. 

In  the  more  general  cases  of  backlash,   the  inertia  and  viscous  fric- 
tion on  the  output  side  of  the  gear  train  are  not  negligible.     For  second 
order  systems,    such  cases  can  be  analyzed  using  phase  plane  techniques. 
This  is  true  for  any  subdivision  of  inertia  and  friction  between  drive 
member  and  driven  member,   for  any  gear  ratio,   and  for  plastic  or 
elastic  impact  at  the  gear  tooth  surface.     While  actual  graphical  integra- 
tion to  construct  the  phase  trajectories  is  too  laborious  to  be  a  tool 
for  generalized  analysis,  understanding  the  techniques,   and  actual  analysis 
of  one  or  two  specific  cases  helps  build  appreciation  for  the  nature  of  the 
problem       Graphical  solution     points  out  the  physical  phenomena  and 
boundary  conditions,   and  leads  naturally  into  the  preparation  of  digital 
computer  programs. 

When  the  backlash  is  taken  up  and  the  system  is  operating  linearly, 
the  equations  of  motion  are  exactly  the  same  as  for  the  no  load  case. 
Considering  here  that  the  motor  used  is  specifically  a  D.    C     armature 
driven  motor,   and  using  a  detailed  notation,   the  equations  of  motion  are: 

—  (  n2  j    +  jt  )  e    +^r  (n2  f    +  £    )  e  =  t  u-n 
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Where  K=  Motor  torque  constant 

=  Motor  generated  voltage  constant 

=  Error  measurement  constant 

=  KlK3/R 

R=armature  resistance 

N=gear  ratio 

J=  inertia  of  system  without  load 

J=    load  inertia 

£=    viscous  friction  (mechanical)  of  system  without  load 

f=    viscous  friction  of  load 

F=    N  f      +  K,   K^  N/R  =  viscous  damping  of  system  without 
m         m  12/  ,  ,     •,.  '         ,  „,,„ 

load  but  including  motor  back  EMF 

FT=f,/N 

bmmanded  displacement 

=  displacement  of  motor 

G=  displacement  of  load 

©=  displacement  of  combined  system. 

To  permit  convenient  numerical  manipulation,   let  N=   1  and  0  =    1, 

R 

then  equation  3-6  becomes: 
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Equation  3-7  may  be  manipulated  into  the  isocline  equation  in  the  usual 

r 

fashion: 
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where   N     is  defined  to  be    d  ©    /d  0     and  is  the  slope  of  the  phase 
1  c/  c 

trajectory  as  it  crosses  the  designated  isocline. 


(3-8) 
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Equation  3-8  is  the  normal  isocline  equation  for  a  linear  second  order 
feedback  control  system.     It  is  expressed  in  terms  of  motor  inertia 
and  friction  and  the  load  inertia  and  friction  because  this  subdivision 
becomes  important  when  backlash  is  considered 

When  there  is  backlash  between  motor  and  load,   and  when  this  back- 
lash is  inside*  the  closed  loop,   the  system  enters  the  backlash  region 
(i.e.  ,   the  gears  lose  contact)  when  the  slope  of  the  load  deceleration 
trajectory  is  exactly  equal  to  that  of  the  linear  system.     In  terms  of 
equations,   when  the  load  is  drifting  free: 

JL®L+FLeL=°  (3"9) 

from  which 

d©  F 

d?r-=    NL=  —it—  (3'10) 

When  the  system  is  operating  linearly  as  described  by  equation  3-8,  and 
the  state  point  reaches  the  isocline  N  which  is  numerically  equal  to  Nt  * 
then  the  gears  lose  contact,  i.e.  ,   the  backlash  separation  line  is  that 

isocline  of  the  linear  system  for  which  N=N    =  —  F       IJ 

1  L  ij  L 

When  the  backlash  separation  line  is  crossed  the  load  drifts  freely 
according  to  the  relationship  of  equation  3-10,  which  guarantees  a  phase 
trajectory  that  is  a  straight  line,     The  motor  shaft  does  not  remain 
stationary.     It  moves,   but  is  driven  open  loop  rather  than  closed  loop, 


*It  is  possible  to  have  backlash  in  the  gear  train  between  motor  and 
load  but  so  located  that  it  is  between  the  measurement  pick  off  and  the 
load.     Then  there  is  no  backlash  in  the  closed  loop,   which  operates 
linearly.     This  case  cannot  have  a  limit  cycle  and  is  not  considered  here 
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with  a  time  varying  input       The  equations  for  the  motor  shaft  motion  is 

j     e     +f     e     =K(&—  e_)  (3-n) 

m      m  m      m  R  L 

The  phase  plane  representation  of  these  equations  is  shown  in  Fig  =    3-1. 
The  separation  dividing  line  is  the  isocline  passing  through  3£.     Equation 
3-11  cannot  be  represented  by  an  isocline  relationship,    so  first  the 
analysis  requires  a  time  solution  of  equations  3-9  and  3-11.     The  back- 
lash is  completely  taken  up  in  the  reverse  direction  when  the  difference 
between  the  displacement  of  the  load  and  the  displacement  of  the  motor 
is  exactlyA,   the  magnitude  of  the  backlash.     In  the  solution  of  equations 
3-9    and  3-11  the  initial  conditions  are  determined  at  the  separation  dividing 
line  and  are  as  shown  on  Fig.    3-2.     Transforming  equation  3-9: 


s     0      (s)   —   s  0  —  © 

L  L  1  L 


+  F 


s  ©        (s)   —  © 
L  J — > . 


=  O      (3-12) 


from  which 


<L 


e 


eL(s)  = 


1 


,y</Ji 


+ 


h 


(3-13) 


Taking  the  inverse  transform 
©_ 


1 


e    (t)  =  e    +  . 

L  Li    fl|jl 


(  i 


Equation  3-11  transforms  to 


FLt 


L         ) 


(3-14) 


m 


s     ©        (s)     —   s  ©         —  © 
m  m,  m, 

1  1 


+  F 


m 


=  K 


s  ©       (s)   —  0 
m  m 


f  -eL(3) 


(3-15) 
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e, 


ISOCLINES  a   TRACT 
ORIES  FOR  THE 
DRIFTING  LOAD 


FIG.  3-1     PHASE    PLANE    REPRESENTATION   OF  SERVO  WITH 

INERTIA   AND   FRICTION  ON  BOTH  SIDES   OF  BACKLASH 
N  =  -f/J=  SEPARATION    ISOCLINE;  A  =  TRAJECTORY 
FOR     COMBINED    SYSTEM;    B= TRAJECTORY     FOR 
DRIFTING    LOAD;    C=  TRAJECTORY     FOR   MOTOR 
AFTER    SEPARATION. 
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TAN¥* 


i-e 


^ 


Ml    . 


I     . 


K 


e   =  e 

Ml         LI 


—  ec-eR-ix> 


e, 


•IG.3-2     RELATIONSHIP     BETWEEN    0LI,eM|10LI    and   0 


Ml 
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since  0      (t)  =  u  (t)  and  ©      (  s)  =   l/s. 
R  R 

Substituting  equations   3-13  for  0      (s),   and  noting  equivalent  initial 

conditons  from  Fig.    3-2,    equation  3-15  may  be  manipulated  to 

0  ©  K  (  1  -  ©         ) 

nrv  m  mi 

°m(s)  =        s         +s(  s  +  F       IJ      )  +  1       s2(s+  F     |  J     ) 

mm'         m  m     m 

K    © 
m 


J       s2  (s  +  FT|JT    )  (s  +  F      I    J     ) 
m  L|     L  n^l      ,m 

Taking  the  inverse  transform  of  equation  3-16  gives: 

J  ©  Ft  K  (1  -    ©         )  J 

mm  m  mm 

0      (t)  =  0         + (   1   -  e~~J  ) — r-, 

m  m  F  m  F      c 

1  m  m 


Ft  Ft  K  0™ 

1 


m  rn  m ,  L 


(1  J  -   e  J 

m  m 


F         F 
L         rn 


(3-17) 


xr  t2  F      t  t2  Ft 

F  FT  J  m  J     T  L 

m  L  m  — L 


(    — 5-     e  J  = e  ) 

J      F  J       F  F2  m  F2  JT 

L      m         m      L  m  L  -L 

If  equations  3-14  and  3-17  are  plotted  with  displacement  as  ordinate 

and  time  as  abscissa  they  may  be  used  as  a  nomograph  to  evaluate  points 

on  the  backlash  dividing  line  for  recontacting  of  the  gear  teeth.     If  both 

are  plotted  on  the  same  coordinate  axes,   then  at  any  instant  the  difference 

between  the  ordinate  values  is  the  relative  displacement     0,-0         , 

L  m 

and  whenever    ©      -     ©        =    Z\    the  backlash  is  taken  up       Since  this  point 
L.  m 

also  defines  values  for      ©        and  for      t,   differentiation  of  the  proper 

equations  permits  evaluation  of        0  and      ©T      at  that  instant. 

m  Li 
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For  repeated  cycles  these  time  curves  need  not  be  replotted,   but  may 
be  used  as  a  nomograph  by  simply  adjusting  the  ordinate  scales.     This 
is  due  to  the  fact  that  the  only  change  in  the  physical  situation  is  a  change 
in  the  initial  conditions,   and  from  Fig.    3-2  there  is  a  fixed  relationship 
such  that 


1   — 0  1-0T 

"1  Ll 

tan»=         -g -     -5 (3-11 

"l  Ll 


where  tan  Wis  known    number.     Then    1   -0         =  0         tan  M*  which    (3-19) 

m  m 

indicates  the  required  adjustment  in  ordinate  scales.     A  pair  of©      and 

©       curves  is  shown  on  Fig.    3-3. 
m 

At  the  instant  backlash  is  taken  up  both  motor  and  load  have  finite 

velocities.     Therefore  both  the  law  of  conservation  of  momentum  and 

the  law  of  conservation  of  energy  must  be  considered.     If  the  impact  is 

defined  to  be  inelastic  (plastic)  this  is  equivalent  to  stating  that  sufficient 

deformation  occurs  to  dissipate  some  energy  as  heat;  then  the  law  of 

conservation  of  energy  need  not  apply  to  the  dynamics  of  the  mechanical 

system,   but  the  law  of  conservation  of  momentum  still  restricts  the 

boundary  conditions.     It  follows  that  after  recontact  of  the  gear  teeth., 

the  requirement  that  momentum  must  be  conserved  is 

J      ©T    +  J         ©       =  (  J       +  J      )  0  (3-20) 

LLmm  m  L,         L. 

where  ©,    is  the  instantaneous  velocity  of  the  recombined  system „     Note 

L2 
that  the  points  at  which  the  gear  teeth  recontact  define  a  dividing  line  on 
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1.5  2  2.5  3 

TIME,  SECONDS 

FIG.  3-31       DISPLACEMENT    CHARACTERISITICS     OF    MOTOR 
AND   OF    LOAD    DURING     BACKLASH    TAKE-UP. 


the  phase  plane,   and  the  result  of  the  applying  the  momentum  restoration 
defines  another  dividing  line  on  the  phase  plane.     These  lines  may  be 
established  point  by  point  in  evaluating  the  phase  trajectory,   but  no    simple 
equation  for  these  lines  has  been  obtained. 

A  phase  plane  plot  for  a  particular  case  is  given  on  Fig.    3-4       This 
plot  shows  the  three  dividing  lines,   and  the  nomograph  of  Fig.    3-3  was 
used  in  computing  Fig.    3-4.     Note  that  the  system  goes  into  a  limit  cycle 
along  path  ABCD. 
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4.     COMPUTATION  BY  DIGITAL  COMPUTER 
4.  1     Introduction 

Use  of  the  phase  plane  techniques  of  section  3  is  instructive  but 
prohibitively  laborious  if  a  large  number  of  cases  is  to  be  investigated. 
Alternatives  are  analog  simulation  or  digital  computation.     The  indivi- 
dual equations  are  readily  simulated  on  the  analog  computer ,  but  instru- 
mentation of  the  boundary  conditions  can  be  difficult.     An  analysis  using 
analog  simulation  has  been  carried  out  elsewhere*  but  was  restricted 
to  the  inelastic  impact  case       It  was  decided  to  use  a  digital  computer 
to  calculate  the  desired  data,   and  the  computer  was  to  be  programmed 
using  phase  plane  concepts  so  that  all  of  the  data  required  for  phase 
plane  plots  would  be  available  when  needed.     The  first  program*  was 
prepared  for  the  inelastic  impact  case  only.     It  operated  successfully 
as  verified  by  point  checks  calculated  long  hand  on  the  phase  plane. 
While  the  results  were  accurate,   only  a  limited  amount  of  data  was 
computed  due  to  lack  of  available  time.     Upon  resumption  of  the  pro- 
ject it  was  decided  to  reprogram  the  computer  to  handle     he  general 
case  of  elastic  *  contact  of  the  gear  teeth,   including   the  case  of 
plastic  impact  as  a  limiting  case       This  program  was  completed 
and  operated  satisfactorily.     It  is  given  in  detail  in  the  appendices. 
The  results  reported  in  the  remainder  of  this  text  is  based  on  com- 
putations made  with  this  program. 


*Knoll  &  Narud 

*New 

*Anderson  &   Luckett 
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TABLE  4=1  SYMBOLS* 


R 


E 
K 


m 


Lm 


P 


m 


m 


A 


N 


(ONE)** 

(THETAL) 
( THE  T  ADL) 

( KM  OT  CONST) 


(RHO) 

(MOINERT) 
(JLOADVAL) 


(FMOTOR) 

(FLOAD) 

(DELTA) 


Angular  step  input  (radians). 
Desired  load  displacement  at  steady 
state 

Output  displacement  of  load  (radians) 

Output  velocity  of  load  ( radians) 

sec 

Error  in  radians.     E=0^   —  © 

R  c 

Constant  of  proportionality  relating 

position  error  to  torque  developed  by 

the  motor . 

Developed  torque  of  motor. 

Torque  of  load. 

Torque  of  load  referred  to  motor. 

Number  of  teeth  on  gear  1 
Number  of  teeth  on  gear  2 

Inertia  of  motor  measured  at  motor. 

Inertia  of  load  measured  at  load. 


Total  inertia    measured  at  motor 
J  =  Jm  +  p2  J 


Friction  of  motor  measured  at  motor. 

Friction  of  load  measured  at  load. 

Backlash  of  gearbox  measured  in  radians 
at  the  output 

Slope  of  phase  trajectory  (combined 
system. 

d©       e 

c  c 


N 


A 


de 


*In  this  table  some  of  the  symbols  are  redefined  and  have  different 
meanings  than  in  preceding  sections. 

*#(  )   Terms  indicate  computer  mnemonics  referred  to  in  Sec.    3, 
Computer  Program  Development. 
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N. 


.   1       .  1 

e   ,  e 

m         c 


f 


w 

£ 


O 

2 

n 


L 


m 


(RESTITUT) 
(ZETA) 
(OMEGAN) 
(OMEGANSQ) 

(FLFTPRIN) 
(INERTRAT) 


Slope  of  phase  trajectory  (load  f lot  ling 
free)  . 

de        e 

N,   * 


L       d0  g 

c  c 

Velocities  of  motor  and  load  after  impaqt 

Coefficient  of  restitution. 

System  damping  coefficient. 

System  natural  frequency. 

System  natural  frequency,    squared. 

Load     riction  ratio     with  respect 
to  friction  of  the  system. 

Inertia  ratio. 


m 


e 


m 


(THETAM) 
(THETADM) 


Displacement  of  motor  (radians) 


Velocity  of  motor  J  radians 

sec 
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4.  2    Development  of  Equations 

The  equations  developed  were  those  of  a  phase  plane  analysis. 
Prior  to  defining  the  net  system  equations ,   the  assumptions  under  which 
the  analysis  was  made  will  be  stated.     When  necessary  these  assump- 
tions will  be  amplified  and  referred  to  later  in  this  work. 
It  was  assumed  that: 

1  .     The  gear  teeth  were  initially  in  contact  and  the  initial 

conditions  of  the  system  were  all  equal  to  zero.     This  was 
later  proved  to  be  an  unnecessary  limitation  for  the  study 
of  steady  state  response. 

2.  Plastic  deformation  of  the  gear  teeth  during  steady  contact 
and  impact  and  any  torsional  deformations  of  driving  shafts 
are  negligible . 

3.  The  inertias  of  the  gears  and  drive  shafts  are  considered 
as  part  of  the  load  or  motor  inertia  depending  on  their 
attachment  in  the  system. 

4.  The  law  of  conservation  of  energy  was  completely  satisfied 
in  only  the  perfect  elastic  case  by  maintaining  the  total 
mechanical  rotational  energy  of  the  system  constant  at  the 
instant  prior  to  and  after  impact.     The  law  of  conservation 

of  momentum  is  satisfied  in  plastic,  elastic,  and  intermediate 
cases.  When  the  law  of  conservation  of  momentum  is  the  only 
equation  required  to  be  satisfied,   the  energy  lost  from  the 
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system  is  dissipated  in  the  heat  of  infinitesimal  deformations 
of  the  gear  teeth 
5.       The  gear  teeth  are  in  contact  only  instantaneously  during 
impact  for  the  elastic  and  intermediate  cases  (excluding 
plastic)  and  that  the  impulse  torques  of  drive,   friction  and 
bearing  supports,    etc.  ,   are  zero  during  impact,   e.g. 

5-     -     ■--■ 


T       dt    =  O  and      G  =  0 

m 


dt-^  o 

6.  The  coefficient  of  restitution  of  the  two  opposing  gear 

teeth  is  the  same  or  is  described  by  an  equivalent  coefficient 
if  the  two  gear  teeth  are  of  unequal  coefficients. 

7.  Backlash  is  assumed  to  be  equal  at  all  points  on  the  gear 
circumference.     Backlash  is  measured  at  the  output  shaft. 

A  block  diagram  of  the  system  considered  is  presented  in  Fig.    4-1, 

and  symbols  have  been  defined  in  Table  4-1  . 

The  equations  for  the  motor  and  load  in  contact  are: 

E  =  ©      —  G  (4-1) 

R  c 

T     =  KE  =  K  (0^  —  0  )  (4-2) 

m  Re 

TT  =  JT    G     +  f_    G  (4-3) 

L         L      c         L,     c 


T 


Lm 


=  pTL  =  p<JLec  +  £Lec>  (4"4) 


0     =  0  G        forward  G     =  rt  G     +   A  backward        (4-5) 

drive  drive 
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Fig.   4-1  Block  diagram  of  second  order  servo- 
mechanism  with  backlash  in  the  gears  between 
motor  and  load. 
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O    =  p  0      forward  and  backward  drive  (4-6) 


m 

hence  T_        =  p  (Jt      0      +  fT    p  ©    )  (4-7) 

Lm       r       Lp      m         L  r      m 

=p2  (JT    ©      +  fT    ©    ) 
r         L     m        L     m 

T       =  J       ©     +  £       ©      +   TT  tA    0. 

m        m      m       m      m  Lm  (4-8) 

j     ©     f     ©       p   (j    ©      f.T  e   ) 

m c_  ,  rn c_     *  L c_      L c_ 

p  p  f  P 

(J  p  JT)  ©  (f  p  f    )  © 

m       r     l'      c        v_m        "    L       c 
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T     =  K  (©      —  ©  )  (4-9) 

m  R  c 

J  £  .  (4-10 

K0=( — +pO      ©     +(—  +pfT    )©+K© 

R         p  r    L  c  p  rL         c  c 


f 

(   —  +  pf_  )  ,  K  © 

e  +lf£ -±—  ©  +  - £~ — =  (4-ii) 

c     (J™  +  PJL)  C       (fm    +PJL) 


KGR 


J  p  J 


m 

+ 


©     +  2fw<3     +  Wn    ©     =U)n    ©D  (4-12) 

r  v        f)     c  C  K 


2 
where  U^n     -  pK 


2  i*Wn  = 


f 


m 


2  (4-14) 

m  '    r     JL 


J        h  p     J 
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The  equation  for  load  alone  with  gears  not  in  contact  is: 

J      0     +  fT    0     =     O  (4-15) 

L      c         L      c 

or  ••      £L     • 

0     +     -rr=—      ©        =    O 
c         J  c 

l_l 

The  equation  for  the  motor  alone  with  the  gears  not  in  contact  is. 

J       G      +f       0      =K(G      —  G)  (4-16) 

m      m        m      m  R  c 

The  method  of  analysis  used  by  New,   Ref.    e,   is  used  in<  order 
to  be  able  to  examine  the  system  independently  of  its  tJ_ 
Defining  first       U)nt  =  t* 


and  differentiating   i^ndt  =  dt* 
dt*  .     i.i  2       .    dt5* 


2 


UJn  =  —  and    ujn    =(-2—) 

dG  dG        dt*         de 

c  c  c 


jfe 


Gc"dt       "     dt*       dt        '   dt*  Gc  (4-17) 

hence      0     =0     U)n  (4-18) 

c         c       " 

similarly 


/dt  *f         d  G 


d   ©  d    G 

*•  c  c 

G     = 


dt 


Idt  f)Z    \dt    /        (dt " ) 


Un  <4"19> 


=  0         and  ©     =  ©       UJn  (4-20) 


(at*) 


2  c  c         c 


finally                  .              #1s  ,  J{ 

"<■•                 ^  at '                                                                      i  a    m 

0  =  ©  U>n=  ©  -rr —                                                       (4-21) 

c         c                c  dt 

JC  •*       *  * 

©     dt  =    \  ©     dt      or  0     =  ©  (4-22) 

J      C  C         C 
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Using  the  equations  of  the  combined  system  and  the  load  alone, 

0  +  2  y OJn0  +    u)n    ©   =  UJ(\    ©  combined  system        (4-23) 

©  +  — - ©    =0  Load  alone  (4-24) 

C       ij    _  c 


and  making  the  indicated  substitution  to  a  transformed  (*)  coordinate 
system  for  the  system  equation 

Wn         &c+z£wn    ©    +Un©c=U)n    ©R  (4-25) 

©c+2fec+e     =eR  (4-26) 

and  introducting  the  slope  equations  of  the  phase  plane 

n   e  +  2  f  e    +  e    =  e„  (4-27) 

s     c         *>      c  c  R 

If  the  load  equation  is  first  put  in  the  form  for  the  phase  plane 

nt  e  +  - —  e   =  o  (4-28) 

L     c      J  c 

^  .*      . 

and  then  transformed  to  the  (*)  coordiante  system  where  0    5t  O 

c 

-*       £L       « 

NT    Q  + 0    =  O  (4-29) 

L     c     J  c 

J_i 

The  same  result  may  he  obtained  by  setting  IjJL^  1  .     The  results  of 
such  a  transformation  require  inverse  scaling  for  practical  applica- 
tion,   examples  of  which  given  later.     It  is  noted  that  the  slope  of  the 
load-free  equation  is  mathematically  the  same  whether  in  the  phase 
plane  or  in  the  transformed  phase  plane.     It  is  pointed  out  at  this 
time  that  the  equations  later  developed  to  satisfy  the  laws  of  con- 
servation of  momentum  and  energy  are  independent  of  the  system 
natural  frequency. 
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To  establish  equations  for  the  law  of  the  conservation  of  momentum 
and  a  relationship  satisfying  the  law  of  the  conservation  of  energy, 
the  schematic     of  Fig     4-2  is  used.     It  is  seen  that 

0   =   —  P^i>    ®   =   —  P^i  ana  since  P   =   1  an  equivalent 

schematic  is  shown  in  Fig.    4-2b.     By  representing  the  inertia  of  the 
load  and  motor  as  inertia  of  the  gears  figure  4-2c  is  obtained  at  the 

instant  of  impact.     The  torque  equation  may  be  expressed 
d© 


t=  j  e  =  J 


dt 


Impulse  is  equal  to  the  time  rate  of  change  of  momentum.     The  ex- 
pressions for  the  rates  of  change  of  momentum  of  the  gears  treated 


separately  may  be  written: 


J       dt-»o 


0 


r       dt    = 


s~       m 
J       m 


d  © 


m 


m 


dt 


dt 


(4-30) 


—  1 


U 


dt       =      J       (0     _  0     ) 
mm  m 


dt-jo 


(4-31) 


i 


and 


F>j 


dt  ->o 


dt      =      J 


m 


(   0       —   ©     ) 

m  m 


for  the  motor 


s 


2J 


r2     dt 


d0 
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Fig.   4- 2a    Schematic  diagram  of  motor  and  load 
inertias  with  interconnecting  gears. 
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Fig.   4- 2b    Schematic  diagram  of  mechanical  system 
equivalent  to  Fig.   4- 2a. 
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n'^ 


Fig.   4- 2c     Gear  teeth  and  loads  involved 
at  instant  of  impact. 
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where  the  primed  velocities  are  those  following  impact. 

If  the  momentum  is  to  be  conserved  in  the  system  the  impulse 

functions  are  equal  and  opposite 

2 


J 


P    J 


m  ,    «  •      .  L 


i 


r(  e    —  e  )  = — r  (  e    —  e>  (4-35) 

1         m  m  '  1  1 


ri  r2 


and 

J 


J2.J-,e"   _©   )   =,©-0)  (4"36> 

J  p  m  m  c  c 


from  substitution  of  the  equations 

ii. 

=   r 
1  2   >       c 


r,    =  r„        0   =     -  p  G 


-    1              1 

e        e, 

m  +        1 

»   i 
(p  G      - 
v          m 

pG       - 

'       m 

»   i 
ec> 

e    +  e, 

m          1 

0 

c 

The  assumption  that  all  other  impulse  functions  are  zero  at  the  instant 
of  impact  is  restated  at  this  point. 

A  definition  of  the  coefficient  of  restitution  will  be  made  for  pur- 
poses of  this  work: 

4 

(4-38) 

It   will  be  shown  that  this  definition  of  e  will  satisfy  the  law  of 

conservation  of  energy.     Consider  the  case  of  e=0, 

O  =  —  (p  G       —  g')   ,     pG      =G  (4-39) 

r      m  c  r  c 

in  which  case  for  plastic  impact  (e=0)  the  gears  are  moving  at  the 
same  velocity  following  impact. 

The  expression  for  the  coefficient  of  restitution  is  examined  for 
the  case  of  perfect  elastic  contact,  e=l,  for  which  the  law  of  conser- 
vation of  energy  is  satisfied. 
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Solution  of  the  restitution  equation  for  e=l  yields. 


.  i  »i 


p  ©     _  e   =  —  p  e    +©  (4-40) 

r      m  ^c  rn  c 

p(e    +  ©'  )  =  e  +  e  (4-4i) 

'         m  m  c  c 

For  conservation  of  energy 

1  .2  «2  . '2  .' 2 

—-(j     e    +j    e  )  =-— (j     ©     +  j    ©  )  (4-42) 

2  mmLc  2mm  L      c  v 


or 

J  2/2  '2  2 

m        •  •  * 

- —    ©     —  e     )  =  (o    -©    )  (4-43) 

JT  m  m  r  c 

and  factoring 

J 

(0     +0     )(©     —0     )  =  (©  +  ©)(©     -0)  (4-44) 

mmmm  cccc 

JL 
The  momentum  equation  which  is  independent  of  e  is  : 

J 

-^     (  ©       -©)=(©-©)  (4-45) 

p  J  mm  c  c 

or  rearranging      J  .,  ...  (4_46) 

0  J  m  m  c  c 

The  factored  energy  equation  is  now  divided  by  the  momentum 

equation  to  yield 


p(©      +©)  =  (©+©)  (4-47) 

•  m  m  c  c 


i 
m  m  c  c 

which  checks  with  the  definition  of  the  coefficient  of  restitution 

for  e=l . 

The  momentum  equation  implied  from  the  impulse  approach  that 

dt -^  o  hence  ©  =  ©      and  ©       =  ©     •     The  definition  of  the  coefficient 
c         c  mm 

of  restitution  is  independent  of  position.     It  may  be  noted  that  the 
equations  for  energy  and  momentum  depend  only  on  the  inertia  dis- 
tribution,  gear  ratio  and  angular  velocities. 
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In  manipulation  of  these  equations  it  is  seen  that  load  inertia 

2 
may  be  transferred  to  the  ©     ,    shaft  by  multiplication  p    ,   the  same 

transfer  may  be  accomplished  with  load  friction.     Since  ©    ,   and  ©    , 

are  related  to  ©     and  ©      it  may  be  reasoned  that  the  results  for 

-Li  J-J 

p=  1  may  be  extrapolated  to  physical  systems  where  p/  1 

Since  a  total  energy  accounting  may  be  made  for  the  system  with 
the  equations  used,   the  rotational  energy  lost  is  attributed  to  the  heat 
of  deformation  of  the  gear  teeth.     This  transfer  of  energy  could  in 
fact  be  determined  from  the  equations  used.     Thus  the  conservation 
of  rotational  energy  e=  1  is  a  special  case  of  a  broad  interpretation 
of  the  law  of  conservation  of  energy. 
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4.3    Computer  Program  Development 

The  physical  equations  of  the  net  systems  were  programmed  for 
solution  using  the  Control  Data  Corporation  1604  high  speed  digital 
computer  utilizing  paper  tape  program  input  and  magnetic  tape  output. 
An  IBM  717  line  printer  was  used  to  extract  data  from  the  magnetic 
tape       The  Control  Data  Corporation  machine  library  and  the  U.   S. 
Naval  Postgraduate  School  computer  subroutine  library  were  used  for 
assembly     Runge-Kutta- Gill  numerical  integration,   and  decimal 
output 

Several  changes  were  made  in  the  forms  of  the  physical  equations 
of  the  net  system  in  order  to  eliminate  duplicate  computing  operations 
and  to  fit  the  equations  to  a  form  suited  to  the  variable  parameters  . 
In  the  table  of  symbols,    computer  mnemonic  (m)  terms  which  are  used 
in  the  assembly  subroutine  have  been  indicated  by  parenthises  and  will 
be  defined  when  encountered. 

The  equation  of  the  system  with  motor  and  load  combined 

gc=  u)n2  eR  -  u)n  gc  -  2  £wn©c  (4"48> 

was  put  in  the  form 


2  ^       2 


G  =U>n  ©„  —  u)^  0  —  A  e        where  (4-49) 


A_  =  2  (^uJ^and 

THETAL  =   B  ONE  =   0 


c 


R 


THETALD  =  U  =0  ZETA  - 
-     c 

UDOT  =  ©  QMEGANSQ  =U)n 
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(RUNGE),   numerical  integration  method       The  RUNGE    sub- 
routine required  the  definition  of  synonymous  terms  for  the  four 
iterative  cycles  used  to  produce  one  extrapolated  set  of  variables. 
The  increment  of  the  independent  time  variable  chosen  was  0.01   sec. 
for  all^O,  1  and  0.004  sec     for  allj*^   0.1. 

The  load  free  equation 

-       £L     • 
was  put  in    he  form      ©    +  r—    0     =  O  (4-50) 

c      J  _         c 
L 

for  the  computer     e    +  B    0     =  O  (4-51) 

c       - —      c 

A  separate  set  of  RUNGE  synonyms  was  used  for  the  load  free  equation 
in  order  not  to  destroy  the  previous  computations  before  they  were 
determined  of  no  further  value  in  computing  and  also  in  the  making  of 
program  decisions- 

©     was  designated  V    and  ©     was  designated  VDOT  for  the  computer 
program . 

The  equation  for  the  motor  alone 

©    +-^  ©      +JL   e  =-|L  e  (4-52) 

m      J  m         J  c        J  R 

m  mm 

was  put  in  the  form 

f 

G      =  ■ -     ©      +    ~      (  0_   -  O   )  (4-53) 

m  J  m        J  Rc 

m  m 

©      =   —  D    ©      +  C    (  ©„   ~  ©  )  (4-54) 

m  —     m      —         R  c 

and  for  RUNGE,   the  terms 

©•      =  THE  T  AM  0     =  THE  T  ADM  =  W 

m  m  — 

and       ©       =  WDOT  were  used. 

m 
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The  equation  representing  the  law  of  conservation  of  momentum 


and  the  definition  of  coefficient  of  restitution  were  combined  to  the 

p  0        (   1  +  e)  +  G     (         ~   ~  et 
rn  c  V        J 


forms         •'  m 

K2       — 


C  T  ~j-"       T 

J-i-pJTL  x  m         '      J 


(4-55) 


m 


i 


G      =0      _epG        +e0 
m      __S V     m         __c  (4.56) 

P 

for  the  computer  program.     No  provision  for  additional  m  terms  was 

made  to  denote  the  primed  values.     The  additional  terms: 

2 
e  =  RESTITUT         n  =    RHO        and  p    =  RHOSQ  were  used, 

Two  equations  were  used  to  define  the  boundaries  of  operation  of 

the  load  and  motor  when  they  were  acting  separatelys 

G    =  n  G        and  G     =    p    G       +       A  (4-57) 

c       r      m  cim  *-* 

Two  major  decisions  of  the  computing  cycle  were  the  determination 
of  the  point  where  the  motor  and  load  would  float  free,  FLOATEST, 
and  the  response  of  the  system  to  impact  of  the  load  and  motor  gears 

when  the  boundary  conditions  on  position  were  met,   COMBTEST  „ 

The  first  decision,   the  point  of  float  free,   is  made  on  the  equality 
or  inequality  of  slopes  in  the  phase  plane ,   i.e.  ,   at  separation 
N    =  N         Equating  the  slopes: 

Li  S  _ 

__LL  /  q,     _=©)—A+B       =0  (4-58) 

©  R  c  — -       — 

c 

The  above  equation  is  solved  after  computing  each  point  in  the 
combined  phase  trajectory.     If    N±  N    J  O  the  system  remains  com- 
bined.    If    Nc  +    N     <  O  the  sytem  separates,   having  the  response 
motor  alone  and  load  alone  immediately  thereafter. 
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The  second  major  decision,   that  of  whether  after  impact  the 
response  should  be  that  of  the  combined  system  or  the  response  of 
the  load  and  motor  acting  separately,   is  made  on  the  basis  of  the 
resultant  velocities  after  impact  and  upon  which  side  of  the  backlash 
the  motor  and  load  positions  are  found. 

The  computer  program,   its  flow  diagram,   a  description  of  the  pro- 
gram operation,   and  an  explanation  of  the  readout  of  information  with 
examples  are  presented  in  the  Appendix. 
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4.  4    Methods  and  Scope  of  Examination 

Two  modes  of  operation  of  the  computer  program  were  used. 
The  first  mode  considered  was  the  phase  trajectory.     This  mode  was 
used  to  check  the  operation  of  the  computing  cycle.     Time,   velocity 
and  position  of  the  load  were  printed  when  the  system  was  combined 
When  the  system  was  separated,   the  values  of  motor  velocity  and  motor 
position  were  also  printed.     Although  the  time  increment  used  for  com- 
putations was  either  0,01   sec.    or  0.004  sec.  ,   printed  outputs  for  the 
first  mode  were  taken  only  every  0.  1   sec.     Additional  printed  outputs 
were  taken  at  the  time  of  contact  of  the  gears  just  prior  to  and 
immediately  following  the  solution  of  the  momentum  and  restitv  In- 
equations.    A  typical  printout  of  the  phase  trajectory  mode  of  operation 
is  shown  in  the  Appendix.. 

The  second  mode  of  operation  was  used  for  pi  ii    ing  only  the  maxi- 
mum computed  positive  load  overshoot  position  for  each  cycle,   the 
associated  load  velocity  and  the  exact  problem  time  of  the  computation. 
A  sufficient  number  of  print-outs  for  each  problem  was  obtained  to 
prove  limit  cycle  existence  and  average  size.     Since  this  investigation 
was  concerned  primarily  with  the  steady  state  response  :,   the  second 
mode  of  operation  was  the  one  utilized  to  obtain  the  majority  of  the  data 
This  read-out  method  markedly  decreased  the  data  reduction  time  for 
the  problem  analysis,   as  opposed  to  the  analysis  of  steady  conditions 
provided  by  a  full  phase  trajectory  print-out       Since  the  time  required 
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for  computer  read- out  was  several  magnitudes  greater  than  the 
computing  time,   valuable  computer  time  was  saved  by  this  mode  of 
operation,     An  example  of  the  printed  output  obtained  when  only  maxi- 
mum load  position  overshoot  was  of  interest  is  given  in  the  Appendix. 
Approximately  1400  phase  trajectories  were  solved  using  this  mode  of 
operation.     Each  solution  required  an  average  of  three  minutes  of  com- 
puter time . 

For  the  purpose  of  general  examination  of  the  problem  it  was 

2 

assumed  that  n  =   1 ,   K  =   1  »U)n=   1,     J     +  p  J    =   1  ,     Specific  solutions 

were  made  with  p  ^  1,  ouL^  1  to  determine  scaling  effects.     11    was 
determined  that  velocities  were  scaled  in  direct  proportion  with  u>n 

and  time  is  scaled  in  an  inverse  ratio  with  u)       It    was  also  determined 

2     fL  2   JL 

that  p      ~ —      andp- —     were  valid  parameters  and  in  agreement  with 

T  m 

the  non  dimensional  parameters  used  for  the  investigation. 

The  program  can  also  be  made  to  solve  linear  systems  by  setting 

A=  o. 

The  major  variable  parameters  used  in  this  investigation  were: 

^=0.1,   0,2,  0.3,  0.4,  0.5,   0.6,   0.8,    1.0 

2f 

— —=  0   ,   0.2S   0.4,   0    6,   0.8.    1  .0 
FT 


2 
L 


P"   J 


=  0.1,   0 .  2  ?   05,   0.8,   0.9 

J 


m      0.9,    0,8,   0.5,   0.2,   0.1 
Z\  =  0  .  3  radians 
e    =  0,   0.6,   08,    10 
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Additional  parameters  of 

p2fL 

— -  =  0.04,    0,1,   0.9;  and^  =  0.01,   0.03,     0.1,0.15  radians  were 
T 
used  in  certain  instances  to  examine  particular  characteristics  of  the 

response . 


The  program  was  limited  somewhat  in  that  the  value  of 
2 
L 


p-J 


=  C^o   and  zero  were  excluded  due  to  generation  of  undefined 
m 
mathematical  quantities. 

The  value  of  backlash  was  made  abnormally  large ,^  =0.3  radians, 

to  allow  an  easier  interpretation  of  the  non- linear  response.     Since 

the  influence  of/\was  linear,   this  caused  no  inaccuracies.     All 

graphs  for  limit  cycle  size  were  plotted  with  ^  =  0.3  radians.     The 

effect  of  values  for  /\  other  than  0,  3  is  also  shown  in  graphical  form. 
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5.      ANALYSIS  CF  RESULTS 

The  primary  purpose  of  this  study  was  the  determination  of  the 
conditions  for  existence  of  a  limit  cycle.     Due  to  the  nature  of  the 
computations  an  analysis  of  the  results  of  each  numeric        case  pro- 
vided an  individual  "yes"  or  "no"  answer,   i.e    ,   a  limit  cycle  existed 
or  did  not  exist.     This  information  plots  conveniently  as  individual 
points,   and  has  been  recorded  on  Figs.    5-1  and  5-2  in  condensed, 
semi -quantitative  form  to  give  a  broad  view  of  limit  cycle  conditions. 
On  Figs     5    1  and  5-2  the  small  circles  merely  indicate  points  at 
which  a  computation  was  made,   the  curved  lines  were  averaged  in 
as  accurately  as  possible  with  the  available  data.     In  general  these 
figures  are  sufficiently  accurate  for  an  engineering  check,    since 
parameter  values  which  are  close  to  a  boundary  obviously  lead  to 
a  long- duration  oscillation  even  if  there  is  no  limit  cycle. 

The  data  in  Figs.    5-1  and  5-2  can  be  presented  in  the  form  of  a 
three  dimensional  picture  which  provides  considerable  perspective. 
Such  a  presentation  is  given  on  Fig,    5-3,  which  is  a  three  dimensional 
sketch  of  the  limit  cycle  boundary  for  the  particular  case  of  e=l  .0. 
The  surface  of  this  contour  is  the  limit  cycle  boundary.     For  com- 
binations of  parameters  specifying  a  point  on  or  below  the  surface 
there  is  no  limit  cycle-     For  combinations  of  parameters  specifying 
a  point  above  the  surface  a  limit  cycle  is  guaranteed,   and  the 
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amplitude  of  the  limit  cycle  is  related  to  the  distance  from  the  surface 
to  the  pertinent  point.     This  is  further  elaborated  on  by  the  sketches 
of  Fig.    5-4. 

For  calculation  purposes  detailed  and  quantitative  curves  are 
desired,   and  the  three  dimensional  figure  is  not  a  convenient  tool. 
Therefore  the  results  of  the  digital  computer  study  are  summarized 
using  a  series  of  families  of  curves.     These  are  given  in  Figs.    5-5 
through  5-28  to  show  the  effects  of  parameter  variation  on  both  the 
existence  and  amplitude  of  the  limit  cycle.     This  is  essentially  the 
data  which  was  summarized  in  Figs.    5-1  and  5-2  (limit  cycle  exist- 
ence) ,   but  is  more  complete  and  more  accurate 

In  all  of  the  computations  reported  to  this  point  the  backlash  was 
set  a^\  =0.3  radians.     It  was  tacitly  assumed  that  all  amplitude  pheno- 
mena would  be  directly  proportional  to  the  amount  of  backlash.     Since 
this  condition  was  easily  checked  by  additional  computations  it  was 

decided  to  verify  the  assumption.     The  results  are  shown  on  Fig.    5-29 

f 
and  5    30,   from  which  it  is  evident  that  the  amplitude  of  the  limit  cycle 

is  indeed  directly  proportional  to  the  amount  of  backlash. 

Included  for  purposes  of  general  information  are  two  phase  tra- 
jectory plots  shown  on  Fig.    5-31  and  5-32.     These  show  the  type  of 
information  available  from  the  digital  computer  data. 
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6.     APPLICATION  OF  RESULTS 

In  this  section,  the  results  are  applied  to  hypothetical  physical 
systems-  Methods  are  indicated  by  which  physical  systems  may  be 
analyzed  or  designed. 

It  is  desirable  to  find  an  equation  relating  certain  system  variables, 
therefore  the  equations; 


2  p    K 


_p_K 


J     +  p  JT        J 
m     r     L         m 


2 

1  +  p    J. 


U 


and 


J 
m 


^wn  = 


£      +  pi  . 
m      r     L 


J       +p   J 

m       ' 


?zi 
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lK 


+   1 
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m 


1  +  p2J 


will  be  manipulated  to  the  form 
f 


m 


£4 


2 
P   f 
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P2f 


+   1 


m         x/ 


P     J 
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or 


2   _    - 


1  +  p     J 


\/p\/K 


m 


A 


m 


+     1 


JmK(1  +  P2jL) 


m 


m 


1/2 


(6-1) 


(6-2) 


—  (6-3) 


(6-4) 


In  view  of  the  form  of  Equation  6-4  for  ^and  recalling  that  the 

2  2 

variable  parameters  used  on  the  figures  were    p  J         and    p   f      = 

L  L 

m  T 
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2 ~    xt  1S  seen  that  the  expression  for     J 

£      +  p    £_  f 

m       '        L  m 

P    fL 

in  fact  includes  the  other  two  variable  parameters. 

The  value  of  overshoot  and  hence  limit  cycle  are  not  affected 
by  Un;  however,   the  time  and  velocity  are  scaled  for  UJ0/  1.0.     Thus 
the  transient  response  characteristics  may  be  related  to  the  same 
parameters  used  in  these  applications . 

If  a  motor  and  load  are  selected  for  a  system  and  the  limitation 
is  imposed  that  system  operation  shall  not  result  in  a  limit  cycle, 
Equation  6-4  and  the  limit  cycle  existence  charts  may  be  used  to 
determine  the  unspecified  system  parameters    3     ,   K  and  p.     The 
coefficient  of  restitution  for  the  proposed  gear  train  material  must 
also  be  determined. 

To  proceed;    select  a  value  of  K  and   p,   then  determine  the  value 

of  y  from  Equation  6-4.     Enter  the  limit  cycle  existence  chart  for 

.  2 

the  value  of  e,    system   y    and  the  parameters    p     J         and 


1 


J 
m 


f 
m 

Determine  if  the  selected  parameters  describe  a  point  of  no  limit 
cycle.     Since  limit  cycle  existence  is  independent  of  backlash  (Zl=  O 
excluded),   backlash  is  not  a  variable  parameter.     Successive  trial 
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values  of  the  unspecified  parameters  may  be  required  to  meet  the 
limitation  of  no  limit  cycle.     If  it  is  aleo  required  that  certain  transient 
response  characteristics  be  met,   these  conditions  could  be  examined 
by  a  similar  method. 

In  the  event  that  no  acceptable  parameters  are  found  to  satisfy 
the  system  requirements  for  no  limit  cycle,   the  figures  for  limit  cycle 
size  may  be  examined  to  obtain  a  minimum  size  limit  cycle.     The 
method  of  solution  is  similar  to  that  previously  described. 

Example:     The  material  intended  for  use  in  the  gears  has  an 

e=  0.6.     The  motor  and  load  have  the  parameters  f     =  0.64,  J     =0.25, 

r  mm 

f     =  0.64,     J    =1.0  and  it  desired  that  the  system  operate  with  no 
L  L 

limit  cycle.     If  the  values  are  selected  K=4.0,     p=  0.5,   the  values 

",      =1.0    and  f        =4.0,      — 
L  m  f 


2                                                                     1 
p  J        =   1.0    and  f        =  4.0, =0.2  result 

T  .  m  f 


m 


m  p     fT  c 

5     3/2 
Solve  Equation  6-4  ior  f  =■    .5  (.64)  (5+1)        ,   thus 

[25(4)  d+li)1/2 
f"  =  0.4.     Enter  the  chart  for  limit  cycle  existence,   Fig.    5-1, 

e=0.6,^=  0.4, =0.2,      p2JL=1.0 

-^    +  1  J~ 

2„  m 

It  can  be  seen  that  limit  cycle  will  exist  for  these  parameters. 

Since  K  is  most  easily  varied  in  equation  6-4s   a  new  K  is  selected 

at  1  .0  in  an  effort  to  find  parameters  for  no  limit  cycle.     It  is  seen 


that   C  -C  \/—    =    0.4    (2) 

JK  =  1.0       J  K=4.0X 
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with  all  other  parameters  constant,   thus$*=0.8         Enter  Fig.    5-1 
with  e=  0  .  6 

£  =  0.8  -±-     =  1.0 


2 

P    J  ! 

- =  0.2 


J  f 

-^+1 

These  parameters  describe  system  operation  with  no  limit  cycle  and 
are  well  within  the  area  limits.     The  transient  response  of  this  system 
may  be  undesirable  and  various  combinations  of    p  and  K  would  have 
to  be  tested  to  satisfy  the  additional  requirements  of  the  problem.     It 
will  be  restated  that  the  existence  areas  of  Figs.   5-1  and  5-2  are  not 
exact. 

The  figures  of  limit  cycle  size  may  be  used  in  the  same  manner  as 

the  existence  charts.     The  same  parameters  of  the  first  case  are 

2                                    1 
assumed  :     k=  4.0,     p  =  0.5,     p  J     =1.0, =  0.2 

J— <  *• 

-— —  m 

J  -r    +  1 

m  2, 

e=0.6,A=0.1   radians.     Enter  Fig.    5- 12  with  the  parameters  given. 

The  ordinate,    Limit  Cycle  Error  (0.3)    =  0.13  radians,   is  obtained. 

A 
The  maximum  positive  error  in  the  limit  cycle  =0.13  (0.  1)  =  0 .  0434 

0.3 
radians.     Enter  the  same  figure  for  e=0.6,^=0.8 

2_  .    «  1 


p  J      =  1,0, 


X  '       f 

m 

J  — +   1 

m  Z. 


=  0.2,  A  =0.5  radians;  this  is  the  second 


L 
case  of  examination  for  existence.     No  y  =  0.8  curve  exists;  hence  no 

limit  cycle  exists  for  this  choice  of  parameters  for  any  nominal  value 
of  backlash. 
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7.      T.F.AN^JPNT  RESPONSE  CHAR/  CTERISTICS;  PT  "  C-TIC  IMPACT 

The  effects  of  backlash  on  the  transient  response  to  a  step  displace- 
ment input  is  also  of  interest,    so  a  relatively  exhaustive  study  was  made 
of  the  plastic  impact  case       The  char    cteristics  of  interest  are  the  peak 
overshoot     Mpt ,   the  settling  time  and  the  peak  time.     Of  these     the  peak 
overshoot,   Mpt,   produced  the  most  satisfactory  pattern  with  the  result 
that  d<  ta  charts  and  equations  have  been  prepared  to  assist  in  analysis 
and  de.;  >gn      Settling  time  and  peak  time  were  much  less  satisfao      y  to 
work  with       The  results  obtained  to  date  are  presented  here,   but  more 
work  remains  to  be  done  . 

The  bulk  of  the  computed  data  assumes  that  the  backlash  is  0.  3 
radians.     To  use  this  data  with  other  values  of  backlash  it  is  first  necessary 
to  establish  the  effect  of  varying  the  amount  of  backlash.     This  is  done  on 
Figs.    7-1,2,3,      These  figures  merely  verify  the  intuitively  obvious 
fact  that  the  amount  of  overshoot  must  vary  linearly  with  the  backlash, 
but  that  the  proportionality  constant  depends  on  the  various  parameter 
values 
7.  1     Maximum  Overshoot;  Plastic  Impact  Case 

Figs.    7-4,  5,6,7,8,   and  9   present  the  results  of  computer  calcu- 
lations in  the  form  of  families  of  curves.     Each  family  is  a  plot   of 

2  2    f     ' 

Mpt  vs  p     J      /  J       for  a  specific  p  /FT     ratio  with  £  as  the  family 

L#     m 


K  &  A  Figs.    5,6,7 

Figs.    10, 11 , 12, 13, 14, 15  in  K  &  A 
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parameter.     All  are  for  the  case  of  plastic  impact       These  curves  can 
be  used  to  predict  the  peak  overshoot  of  any  arbitary  second  order  system 
with  backlash  and  plastic  impact  by  reading  off  the  pertinent  values  and 
inserting  in  the  equation 


Mpt 


A=A1 

^l 


v  A=0 


A 


0.  3 


Si 


A--.3 
S-   Sj 


—    Mpt 


A=0 


(7-1) 


As  was  pointed  out  earlier,   the  solutions  were  made  and  the  results 
were  plotted  using  a  backlash  angle  of  0    3  radian.     Thus  the  valaes  of 
amplitude  of  maximum  overshoot  on  the  curves  are  significant  only  in  that 
they  represent  the  slope  of  the  overshoot  versus  backlash  straight  lines. 
This   slope  must  be  calculated  in  order  to  scale  the  value  of  overshoot  read 
from  the  curve  for  the  magnitude  of  backlash  actually  expected       Therefore, 
in  order  to  refine  the  prediction  scheme,    other  design  curves,   figures  7-10 

through  7-15,   were  calculated  using  the  curves  of  figures  7-4  through  7-9 

dM 
just  discussed.     Along  the  curves  the  slope        pt   was  calculated  by, 

dM  r~x_ 

pt  M  M 

Pt,  -         Pt, 

?  =  ?1  ?=£l  f=fl 


dA 


(7-2) 


These  values  of  dM      /dA  were  then   plotted  versus  inertia  ratios  at 

pt 

constant  friction  ratios  as  before.      These  curves  give  a  more  straight- 
forward prediction  method  with  the  following  equation- 
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=  M 


A=A, 


Pt 


+    dM    J  d/\ 
pt 


M\ 


x     A 


(7-3) 


In  order  to  demonstarte  the  procedure  for  using  the  design  charts 
shown  in  figures  7    4  through  7-15,   an  example  problem  is  solved  obtaining  the 
maximum  overshoot,   M      ,   for  a  second  order  servo  with  backlash.      The 
differential  equation  for  the  combined  motion  of  an  armature  controlled 


motor  is: 


e   + 

c 


f   +  f 
t      1 


J  +p  J 

m  l 


e   + 

c 


PKa 


J       +f>  J 


e   = 

c 


PKa 


(7-4) 


m 


J     +PJ 
L  m    r      L 

Given  the  following  specifications  of  a  D.   C,    shunt  motor  coupled 

through  a  gear  box  to  a  load: 

Horsepower      =      .008  ( 

Rated  RPM        =     4000 


R 

K 

g 

J 
m 

f 
m 

K 
e 

P 

A 

j 

L 

£L 
K 


282  ohms 

=     .  233  volts/ radian/ second 

-6  2 

=     8.2  X  10        slug   feet 

-     3.34X10"      lb. -feet/ radian/ sec. 
=     10  volts/ radian 
=    1 

=       05  radians     K 

-6  2 

=     16.4  X  10        slug-feet 

=     10.4X10"      lb. -feet/ radian/ sec. 

=     .  1464  lb. -feet/ampere 
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The  required  values   shown  in  Equation  7-4    are  obtained  as  follows- 

K  _  t       e  .  1464  x  10     _    .0052    LB-FT 

a  R  282 


/JDK 


U)  =    H-± ./(I)     -0052  ) 14   5       M5 

tl    /  2  /  A  5 a        SEC 

Y  J    +p     J         V8.2xl0"     +(l)M6.4xlO"b 
m  '  L 

f<-=PfT+£       =  t1)      10.4  x  10"     +3.34xl0"6 

6  LB-FT 


13.74  x  10 


RAD/SEC 

,    LB-FT 


FT=£T  +  £lS  (13.74+  121)    x  10~6  =  135xl°         RAD/SEC 

K    K  _6        LB-FT 

f  t      g         _        1464  x  .233  121  x  10  RAD/ SEC 

1  R  282 

Jm+fJL      V     2UV 

2  6 

2  _     .   _  (1)      16.4  x  10"  2.0 

PJl/J-=_8.2x10-6  = 

P  __L  10.4  x  10  "  .077 

FT  135    x  10  "6 

Two  methods  are  shown  for  obtaining  the  maximum  overshoot  with  a 
backlash  angle  of   .05  radians  in  the  gear  trains. 

The  first  method  is  accomplished  using  the  curves  of  figures  7-4 
through  7-9       The  procedure  is  simple,   but  requires  more  steps  than  does 
the  second.     Also  one  must  remember  that  these  curves  are  drawn  for  a 
backlash  of  0    3  radian  and  must  be  scaled  for  the  expected  amount  of  back- 
lash.    The  second  method  makes  use  of  the  design  charts  of  figures  7-10 
through  7-15. 
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Method  1 

Entering  figure  7-4  (p     £      /F      =  0  .  1)  ,   at    />      JT  /  J     =2.0,   and  a 

L  /       r  (  L      m 

damping  ratio  ,y  =  0.19,   the  amount  of  overshoot  for  a  backlash  angle  ,  Z-^  = 


0.3  rad.   is  M       =1.658.     Then: 
pt 


M  =   Mnt 

pt  Pl 

A=-05 
f=.19 


+    .05/. 3 


^=0 

f  =    19 


From  second  order  linear  characteristics 


M 


M 


pt  pt 

A=3 
£-.19 


4=0 

{•=.19 


M 


Pt 


1.545 


hence , 


M       =   1 .545    +   .05/ . 3        1.658      -1.545 


pt 


=  1 .545  +     .0189 


=  1  .564 

Method  2 

M 
pt 

^=.05 
f  =  .19 

=  M 
pt 

+  dM     /dA 
A  =  0              *' 

r=.i9 

x  A 


p.  19 


2  2 

from  figure  7    10,(p     £L/   ^T  =  0.1),   at  JT/Jwi  ==  2.0,   and  a  damping 


L      m 


ratio      j"*=0.19,       dM     /dA    =0.357 

pt 


Also: 


M 


pt 


1.545 


A=0 
f  =.19 


hence , 


M        =  1  .545  +   .357  X  .05 
pt 

=   1 .545  +   .018 


1  .563 
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The  procedure  for  using  the  charts  remains  unchanged  for  servo 
systems  which  utilize  other  types  of  motors.     The  only  difference  which 
arises  is  in  the  calculation  of  the  chart  parameters  from  the  known 
system  constants.     For  a  field  control  motor  if  the  field  time  constant 
is  negligible,   the  equation  is   similar  to  the  armature  control  motor  equation 

J    +p  JT  J    +pj,  J   +p  JT 

m   *        L  m     I       L  m   r      L 

In  like  manner  the  equation  for  the  system  with  a  two  phase  motor  is 

F        K  aK  pK 

e     + 0     +    — — -      0     = - ©  (7    6) 

C  T  <i  c  _      ,    <iT  c        -      .     ^   -  K 


J     +  p  J  T  J     +DJ7  J+pJ 

m    '        L  m  T      L  mTL 


where 


AX 

K     =     -    -r (from  characteristic  curves  of  two  phase  motor) 

n  on 


s)  T 
K       =      — ^ (from  characteristic      curves  of  two  phase  motor) 

1 


de 


The  calculation  of  the  chart  parameters  is  a  simple  matter  once  the  defin- 
ing equations  are  known. 
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7  .  2    Maxim  urn  Oye  rsho  ot ,   E  la  stic  Impact  Case. 

Investigation  of  a  limited  number  of  cases  reveal  that  the  amount  of 
overshoot  with  plastic  impact  is  always  greater  than  that  with  elastic 
impact.     The  difference  between  this  maximum  overshoot  decreases  as 

X 

increases.      However,   under  no  circumstances  was  the  effect 


FT 


onmaximum  overshoot  greater  than     .01   radians.     The  effect  on  the  time 
of  maximum  overshoot  appears  negligible.     In  a  few  cases  the  time  was 
different  for  the  plastic  and  the  elastic  cases  but  the  difference  was  very- 
small  and  no  pattern  of  difference  could  be  ascertained. 

These  results  point  out  the  validity  of  the  original  assumption  made 
for  this  study.     In  a  practical  system  there  will  exist,    of  course,    some 
amount  of  gear  bounce  when  the  gears  contact  each  other.     The  effect  on 
maximum  overshoot,   however,   is  so  small  that  the  curves  obtained  in 
this  study  remain  reasonably  accurate 
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7  .  3    Settling  Time 

In  determining  the  effect  of  backlash  on  the  settling  time  of  a  second 
order  servo  system,   it  was  necessary  to  first  find  a  definition  for  settling 
time  which  would  be  suitable  for  all  of  the  cases  which  were  investigated. 
Normally,    settling  time  is  defined  as  the  time  for  the  output  to  reach  98% 
of  it  final  value.      This  definition  has  real  meaning,   however,    only  for 
linear  stable  systems  in  which  there  are  no  sustained  oscillations.     For 
systems  which  do  exhibit  small  sustained  oscillations,   the  settling  time  is 
sometimes  defined  as  the  time  required  for  the  amplitude  of  oscillations  to 
decrease  to  within  2%  of  the  final  value.     For  very  small  amplitudes  of 
oscillations  this  definition  gives  an  unrealistic  value  for  settling  time. 
For  the  purpose  of  this  study  settling  time  is  considered  as  the  time  re- 
quired for  the  system  output  to  settle  to  98%  of  the  difference  between  the 
magnitude  of  the  input  and  the  amplitude  of  the  sustained  oscillations. 

For  all  cases    for  which  the  damping  coefficient  was  less  than  0.6 
there  were  no  conclusions  as  to  effects  of  backlash  on  the  system  settling 
time  since  the  magnitudes  of  the  sustained  oscillations  were  always  quite 
large  for  the  amount  of  backlash  used  in  the  solutions  and  the  settling  time 
was  difficult  to  obtain  accurately.     For  those  cases  for  which  the  damping 
coefficient  was  greater  than  0.6,   the  sustained  oscillations  were  small 
if  they  existed  and  an  analysis  revealed  the  following  general  conclusions: 
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1  .     As  the  size  of  the  backlash  angle  increases,   the  settling  time 
will  increase  provided  the  other  parameters  are  kept  constant.     For  the 
cases  investigated,   this  increase  in  settling  time  was  nearly  proportional 

to  the  backlash  angle. 

2 

2  As  the  ratio  o   J,  /J       increases,   with  other  parameters  constant 

*        L     m  r 

the  settling  time  increases. 

2 
3.     At  the  ratio  p  f      /   F        is  increaseo,   with  other  parameters 

constant,   the  settling  time  decreases. 

The  data  included  in  Tables  I- III  show  these  general  conclusions 

and  also  some  quantitative  effects  of  backlash  on  settling  time. 
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TABLE  I 
SETTLING  TIME  DATA 


/ 

TO  _ 

Radians 

Seconds 

.6 

1/9   . 

.2 

.3 

28 

.6 

1/9 

.4 

.3 

21 

06 

1/9 

.6 

.3 

16 

.6 

1/9 

.8 

.3 

13 

.6 

1/9 

1.0 

.3  .--■■ 

7 

.6 

1/4 

.2 

.3 

28 

.6 

iA 

.4 

.3 

24 

.6 

1/4 

.6 

.3 

16 

.6 

1/4 

.8 

.3 

14 

.6 

1/4 

1.0 

.3 

9 

.6 

1 

.2 

.3 

29 

.6 

1 

o4 

.3' 

28 

.6 

1 

.6 

.3 

23 

.6 

1 

.8 

.3 

17 

.6 

1 

1.0 

.3 

10 

.6 

4 

.2 

.3 

29 

.6 

4 

.4 

.3 

28 

.6 

4 

.6 

.3 

26 

.6 

4 

.8 

.3 

23 

.6 

4 

1.0 

.3 

10 

.6 

9 

.2  . 

.3 

29 

.6 

9 

.4 

.3 

29 

.6 

9 

.6 

.3 

26 

.6 

9 

.8 

.3 

23 

.6 

9 

1.0 

.3 

10 

110 


TABLE  II 


SETTLING  TIME  DATA 


LJ 


n 


1.0 


f 

1% 

M 

FT 

Radians 

Seconds 

•6 

1/9 

.8 

.03 

5 

.6 

1/9 

.8 

.10 

7 

.6 

1/9 

.8 

.15 

8 

.6 

1/9 

.8 

.30 

13 

.6 

1 

.8 

.03 

11 

.6 

JL 

08 

.10 

14 

.6 

,8 

.15 

16 

.6 

1 

.8 

.30 

17 

.6 

4 

.8 

..03 

12 

.6 

4 

.8 

.10 

14 

.6 

4 

.8 

.15 

17 

06 

4 

.8 

.30 

23 

.6 

9 

.8 

.03 

13 

.6 

9 

.8 

.10 

16 

.6 

9 

.8 

.15 

18 

.6 

9 

.8 

.30 

23 
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TABLE  III 

SETTLING  TIME  DATA 

LJr>  =K0 

• 

y 

/jM 

FT 

Radians 

Seconds 
t 

.8 

1/9 

.2 

.3 

8 

.8 

1/9 

.4 

.3 

7 

.3  . 

1/9 

.6 

.3 

6 

.8 

1/9 

.8 

.3 

5 

.8 

1/9 

1.0 

.3  " 

5 

.8 

•  1/4 

.2 

.3 

20 

.8 

1/4 

.4 

.3 

9 

.8 

1/4 

.6 

.3 

7 

.8 

1/4 

.8 

.3 

6 

.8 

1/4 

1.0 

.3 

5 

.8 

.8 

1 
1 

.2   . 
.4 

.3 
.3 

43 
31 

;  '.<j ;, 

.8 

1 

.6 

.3 

10 

.8 

1 

.8 

.3 

8 

.8 

1 

1.0 

.3 

6 

.8 

4 

.2 

.3 

48 

.8 

4 

.4 

.3 

36 

.8 

4 

.6 

•3. 

12 

.8 

4 

.8 

.3 

9 

.8 

4 

1.0 

.3 

6 

•8 

9 

.2 

.3 

52 

.8 

9 

.4 

.3 

43 

.8  . 

9 

.6 

.3 

13 

08 

9 

.8 

.3 

9 

.8 

9 

1.0 

.3 

6 
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7.4    Time  of  Maximum  Overshoot 

For  each  of  the  cases  investigated  the  time  of  maximum  overshoot 
(peak  time)  was  recorded  from  the  printout  of  the  solution.     It  was 
found  that  the  presence  of  backlash  causes  a  slight  increase  in  the  peak 
time  over  the  linear  system  for  those  cases  where     J    "^.5,   and   has 
negligible  effect  in  those  where     J)    <^.  5  .     When  p  f       is  increased  from 

zero  to  0    5     other  parameters  constant  andj    "^  5  peak  time  increases 


slightly  and    remains  essentially  constant  for    O   f        ^5.     As 


2. 

F' 


PJ     /  J       is  increased,   other  parameters  constant  and  (/">.  5  ,   peak  tim 
L      m  -* 


is  decreased 
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7.5  Other  System  Characteristics  . 

A  careful  study  of  the  results  as  plotted  in  figures  7-4  through  7-9 
reveals  some  interesting  phenomena.     At  all  friction  ratios  it  is  observed 
that  as  the  ratio  of  load  inertia  to  motor  inertia  is  made  smaller  the  ampli- 
tude of  the  overshoot  is  less.     There  is  a  point  where  the  overshoot  is  identi- 
cal to  that  of  a  linear  second  order  system,   in  other  words,   where  backlash 
has  no  effect  on  the  overshoot.     Note  that  this  point  occurs  at  different 
inertia  ratios  for  each  value  of  damping  ratio.     Also  this  point  varies  de- 
pending on  the  value  of  friction  ratio       If  the  inertia  ratio  is  decreased  beyond 
this  point  then  the  amplitude  of  overshoot  remains  constant  at  the  value  of 
the  linear  case.     Figure  7 -  1 6  is  a  plot  of  the  boundaries  of  the  combinations 
of  variables  which  will  cause  an  effect  on  the  amplitude  of  overshoot  if  the 
system  has  any  appreciable  amount  of  backlash       The  coordinates  are  fric- 
tion ratio  and  inertia  ratio  with  the  boundaries  labeled  with  values  of  damping 
ratio.     Areas  to  the  right  of  a  boundary  are  combinations  of  the  variables 
which  affect  overshoot.     These  curves  were  plotted  using  information  deducted 
from  figures  7-4  through  7-9-     Three  points  were  selected  as  test  points  to 
prove  the  validity  of  the  boundaries  as  plotted.      These  points  are  shown  on 
the  figures  as  points  "A"  ,   "B",   and  "C".     The  system  parameters  used  at 


Point  B  Point  C 

f     =0    6  V  =0.3 

2.     ,„  2, 


the 

test  points 
Point  A 

:  were 

T 

.8 
8 

PfL/FT=0.6  P  fL         =0.4 


FT 


P2jL/Jm=i.o  p2VJm=0-25         p2w°-ui 
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The  system  was  solved  on  the  digital  computer  to  obtain  the  maximum 
overshoot  for  these  three  sets  of  parameters.     According  to  the  curves 
point      A    falls  to  the  left  of  the  boundary  of  the  $*  =  .8  area.      Therefore, 
the  overshoe  t  should  be  the  same  as  that  of  the  linear  case.     The  com- 
puter solution  gave  an  overshoot  of  1.019,   which  is,   in  fact,   identical 
to  that  of  the  linear  case.     Similarly,   point    B    falls  to  the  right  of  the 
y    =0.6  boundary  and  in  the  area  where  the  amplitude  of  overshoot    should 
not  be  the  same  as  for  the  linear  case.     The  computer  solution  gave  a  maxi- 
mum overshoot  of  1.110  which  is  slightly  greater  than  the  linear  case. 
Point    C    falls  to  the  right  of  the  J    =0.  3  boundary  and  m  the  area  where 
the  amplitude  of  overshoot  should  be  different  from  that  of  the  linear 
system.     The  computer  solution  gave  a  maximum  overshoot  of  1  .  41  while 
that  of  the  linear  system  is   1  .372.     From  these  boundary  curves  one  can- 
not determine  the  exact  effect  of  backlash  on  the  maximum  overshoot  of  a 
system,   but  can  determine  in  which  systems  backlash  will  cause  an  appreci- 
able effect, 

Another  phenomena  which  is  revealed  in  the  results  as  plotted  in 
figures  7-4  through  7-9  is  that  for  all  cases  the  greatest  maximum  over- 
shoot occurs  in  the  range  0  .  1  (  p   J     /  J       /    1 0  .  0    for  the  range     of 

2 

rj     /J       tested       However,   it  appears  that  since  the  curves  are  trending 
L      m 

2 

upwards  at  the  high  end  of  the     p  J     /J         axis  that  at  very  large  inertia 

ratios  the  overshoot  would  be  greater.     It  is  felt  that  these  large  inertia 
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ratios  are  outside  the  practical  range.     An  attempt  was  made  to  explain 
these  phenomena  but  without  success.     It  is  often  very  difficult,    if  not 
impossible,   to  find  a  physical  explanation  for  certain  behavior  caused  by 
non-linearities        They  can  merely  be  pointed  out  from  analyses  of  experi 
mental  information  and  accounted  for  accordingly 

It  is  interesting  to  note  that  when  the  ratio  of  load  friction  to  motor 
friction  is  below  0.4,   as  in  figures  7-4  and  7-5  and  1 1  ,    there  is  a  possi- 
bility of  an  overshoot  for  the  critically  damped  case,   i.e.    3      =1.0, 
depending  on  the  inertia  ratio.     At  first  glance  one  would  reason  that  if 
the  system  is  critically  damped  there  should  be  no  overshoot.     However, 
it  can  be  seen  that  if  the  friction  ratio  is  small  and  the  inertia  ratio  is 

large  then  the  ratio  of    F     /J     ,   which  is  the  slope  of  the  drifting  load 

J-i       L 

trajectory,   is  made  small.     Therefore,   when  the  load  separates  from  the 
motor,   its  velocity  remains  fairly  constant  at  a  relatively  high  value  and 
its  momentum  carries  it  past  the  correspondence  point. 
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8.     CONCLUSIONS 

Since  this  report  is  to  be  the  first  of  a  series,   it  is  apparent  that 
the  results  are  incomplete  and  the  conclusions  must  therefore  be  brief: 
1  .     It  has  been  established  by  calculation  that  a  second  order  feedback 
control  system  with  backlash  may  or  may  not  exhibit  a  limit  cycle. 

2.  The  existence  of  a  limit  cycle  as  a  function  of  various  parameters 
has  been  established,  and  a  series  of  charts  suitable  for  analysis  and 
design  have  been  prepared. 

3.  The  amplitude  of  the  limit  cycle  depends  on  the  amount  of  backlash, 
and  charts  illustrating  this  dependence  have  been  prepared. 

4.  Peak  overshoot  depends  on  the  system  parameters  and  the  amount  of 
backlash.     Charts  and  eque  tions  suitable  for  analysis  and  design  have 
been  prepared 

5.  The  effect  of  elastic  impact  at  the  gear  teeth  is  si;: lificant,   but  much 
less  so  than  was  anticipated.     Most  of  the  design  charts  include  this 
effect. 
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APPENDIX    A 

COMPUTAT/O/V     FLOW  O/AGGAM 
CDC   /604      D/G/rAL     COMPUTER 


FAULT 


TART 


aoAWONo] ' 


t»  srornow:  88 


t»co*re3T.39 


7*7.7/ 
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from  MOTOfK.38 


COMBTeSTTSA- 


jcOMPLMMr 


55 


57 


56 


% 


I 

5TOP 


pRgPlANS, 


PLU3P0KC 


68 


Yea 


Is. 


4? 


TS 


',67 


t«» 


Xt5 

PL 

1" 

■ 

-»SCO*.R. 

' 

70 

71 

1 

U  LOAO.S7 

t» 

COM  aitjeo. 


fr*m7*,7S.T6 
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front  S9,  SO 


from  27,43 


MINUTMTJhte 


umitcvc 


YES 
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APPENDIX  A  (continued) 
EXPLANATION  OF  FLOW  DIAGFAM  BLOCK  SYMBOLS 

.1.  Clear  fault  stop. 

2.  Print  last  overshoot  stored  in  LBGTBUF  (symbol  888). 
Print  out  position  in  phase  trajectory  of  fault  stop 
(symbol  666) . 

3.  End  of  file. 

k.     Actuate  fault  stop. 

5.  Load  zeta  from  Index  1  (b  ). 

6.  Print  last  overshoot  stored  on  LTMTTBUF  (symbol  888). 
Print  position  in  phase  trajectory  where  it  was  manually 
stopped  (symbol  lft6l). 

7.  Stop  if  Zeta  =  0. 

8.  Zf   0J„  =  a 

9.  Is  zeta  =  0.1? 

10.  Is  zeta  less  than  0.1? 

11.  Set  up  computation  time  of  0.01  seconds  for  problem  and 
set  up  for  print  at  every  0.1  seconds  when  zeta  is  greater 
than  0.1. 

12.  Set  up  computation  time  of  0.00*1-  seconds  for  problem  and  set 
up  for  print  at  every  0.1  seconds  when  zeta  is  equal  to  or 
less  than  0.1. 

13.  Set  up  limit  cycle  print-out  routine  for  20  overshoots  and 
average  the  last  8. 

l^.  Set  up  number  of  limit  cycle  print-outs  and  averaging  routine 
depending  upon  the  size  of  zeta: 
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f 

Overshoots 

Start 
Averaging 

ULvide 
by  * 

Comp. 
Time 

■  ^0.8 

6 

1* 

2.0 

0.01 

0.6 

T 

5. 

2.0 

0.01 

0.9 

8 

9 

3.0 

0.01 

O.k 

10 

7 

3.0 

0.01 

0.3 

lif 

10 

if.O 

0.01 

0.2 

20 

12 

8.0 

0.01 

^0.1 

17 

12 

5-0 

0.00*4- 

*The  last  row  printed  has  this  number  printed  on  the  far  left 
if  the  problem  has  a  normal  solution.  This  number  indicates 
how  many  overshoots  were  averaged.  The  average  value  has  1.0 
subtracted  from  it  and  is  printed  in  the  far  right  of  the  same 
last  row. 

15.  Page  reject. 

16.  Clear  buffers  and  set  up  for  new  problem. 

17.  Restitution,  e,  from  Index  3  00. 

3 

Backlash,  ^  ,  from  Index  \   00. 

Jm  from  Index  5  (b_ ) . 

5 

*    .     Jt  -Jm         K     „     Jrtx 

*  Jl>  — : '  m   C 


/*■  **. 


j 


&-    from  Index  6  (bg). 

*       _/»**  *  ft,    -    *j£  -  D 

18.  .  Have  all  values  been  cycled  through? 

19.  Set  ST0FN0W  after  this  run  to  stop  absolutely. 
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20.  Backlash,  ^  ,    =0? 

21.  Print  consecutive  number  of  this  parameter  run  and  all  the 
parameters  used: 


*/>«*•  A         J>0t 

Vpt 

e           A 

/* 

Jm              Jl 

JmV°£Ji. 

fm           fu 

/>* 

K        KAr  c 

*fco 

VjjB    6>* 

jr 

22.  Set  program  so  always  acts  as  combined  system. 

23.  Print  consecutive  number  of  this  parameter  run  and  selected 
parameters  for  linear  systent 

2k,     Solve  the  second  order  equation  by  Runge-Kutta-Gill  numerical 

integration  every  0.01  seconds  problem  time  when  zeta  is  greater 
than  0.1,  otherwise  every  0.00!+  seconds.  Initial  conditions 
are  ©0  =  1.0,   ©  =  Q-   "  0  said   (^  m  1.0  or  the  pre- 
viously  computed  point.  Next  point  is  computed  by  four 
iterations  of: 

ec=  cott(©R-0c)-&/coh  ©c 

25.  Set  EXIT  from  print  routine. 

26.  Load  time,   ^c  and  Sfc  in  Print  "buffer. 

27.  Go  to  LIMITCYC  to  print  only  the  overshoot.  Otherwise  print 
the  phase  trajectory  points  for  each  0.1  seconds  up  to  a  maxi- 
mum of  12  cycles. 

28.  Store  0  for  Q^  and  ©m 
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29.  Is  problem  time  0.1  sec? 

30.  Is  problem  time  greater  than  0.1  sec? 

31.  Print  the  number  of  point  computed,  time,  ©c  and  (3C  . 

32.  EXIT  (automatically  set  for  desired  jump-out. ) 

33.  Store  time,  Oc  7    ®c   and-  ©c  for  mator  an<^  i-08-*1  initial 
conditions. 

?k.  Is  slope  1^  m   N  ? 

-r-* —   ^Jr^ta"T  -!-=-  m  positive  value? 

35»  Is  (^  negative? 

36.  Store  0C-  A     in  ©  m 

37.  Solve         <3C=    -Vjj^9c)   f°r    ©c         and    &C 

by  Runge-Kutta-Gill. 

for      ©*>      and      ©^by  Runge-Kutta-Gill 

39.  ©L  -/*  em  ^  o  ? 

40.  0  l   ""Z0  ©  h*  ""  A     =  negative? 

41.  Set  EXIT  from  print  routine. 
k2.  Set  EXIT  from  print  routine. 
43.  Load  time,  @    and  Qc    , 
kh.  Load  ©M  and  @^  . 

k-5.     Print  number  of  point  computed  and  stored  values. 
k6.     EXIT  (automatically  set  for  desired  Jump-out). 
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V7.  Solve 


g  _.   ____      L©m(.  +  e)+6cOVi-e) 


Y 

q^  ______   ________ 

g'    «-"-6_.&+      K     (Sk-  ©c) 

i*-8.  Set  EXIT  from  print  routine. 

o 

k9.  Is  Qc   positive?' 

50.  Set  up  to  count  one  phase  trajectory  when  ®     is  negative. 

51.  Is  Oc  >  1.0? 

52.  Have  12  cycles  of  phase  trajectory  "been  completed? 
53  •  Arrange  for  no  counting  until  ^c  is  positive  again. 
5k.  Is  /?®y*\"~   ©c   negative? 

55.  Set  up   —  (/°<S>*\  "  ^c)  . 

56.  Is  v00m  -  (_?£   greater  than  preselected  £  ? 
57*  Set  up  combined  system. 

58.  Is  0C  positive? 

59.  Is  Q     negative? 

60.  Is  _PC—  J  positive?  ' 

61.  Is  ^c  —  I  negative? 

62.  Is  <_)c*~  I  negative? 

63.  ©c  -s*Q^  0? 

#*•  ©c"y*^lna*  0? 

65.  Is  <_=»c  -  so  &*i   positive? 

66.  Is  9C  -  yO  6?w»Positive? 
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67.  Is  Slope  NL  =  Ng  ? 


.■■, 


z  H 


J         h  +         m  positive  value? 


68.  Modify  FLQATEST  to  keep  system  combined.  Takes  off  modification 
when   ©e   is  negative. 

69.  Is  Qc  -  /9  ©^  positive? 

70.  Store  €>u-  A   in  ©  ^ 

•° 

71.  Store  Qc   in  ©^ 

»72.     Is    0C   »    I     negative? 
75.    ©c  -•*©**  =  0? 
7^.    Oc"  /*^w  =  0? 

75.  Is  Gc   -/0  ©^positive? 

76.  Is  @c  -yO  ©^positive? 

77.  Modify  FLQATEST  to  keep  system  combined.  Takes  off  modification 

c 
when  ©c  ^s  positive. 

78.  Have  1  l/2  minutes  of  real  time  elapsed  since  the  last  print-out? 

79.  Is  0c  <  l-0? 

80.  Set  up  to  print  vhen  in  the  fourth  quadrant. 

81.  Check  if  in  the  fourth  quadrant. 

82.  Is  the  new  ©c  greater  than  the  previous  ©^  ? 

83.  Store  the  new  time,  ©^  an^  ®C   * 

8U-.  EXIT  (automatically  set  for  desired  Jump-out),. 

85.  Print  the  last  overshoot  stored  in  LBffTBUF  (symbol  888) 

86.  Stop  real  time  clock. 

87.  Print  position  in  phase  trajectory  where  it  was  stopped  by  the 
real  time  clock  (symbol  999). 
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88.  Return  to  "beginning  to  start  new  problem  if  new  parameters 
remain  to  "be  evaluated.  Unconditional  stop  if  all  parameters 
nave  been  evaluated. 

* 

89.  Has  the  required  number  of  overshoots  been  evaluated  prior  to 
starting  the  averaging  routine? 

90.  Add  the  ^c  just  evaluated  to  the  previous  sum  for  averaging. 

91.  Add  one  to  the  count  prior  to  averaging. 

92.  Clear  time  clock. 

o 

93.  Print  time,    ©c    and  O^  . 

$k.  Required  number  of  total  print-outs? 

95.  Clear  buffers. 

96.  Obtain  the  average  of  G£  . 

97.  Subtract  1.0. 

98.  Print  the  average  of  ($^—  1.0  and  how  many  overshoots  were 
•  utilized  to  obtain  the  average. 

99.  Is.  ©c  less  than  1.005  radians? 

100.  Clear  buffers. 

101.  Start  real  time  clock. 

102.  EXTT  (automatically  set  to  desired  jump-out). 

103.  Has  30  seconds  of  real  time  elapsed  since  the  last  print-out? 
lO^f.  Print  the  position  in  the  phase  trajectory  where  it  was  stopped 

l>y  the  real  time  clock  and  Ql   being  less  than  1.005  radians 
(symbol  1003 ) • 
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Jump  Switches: 


J  -  For  LIMITCYC,  allows  the  maximum  overshoot  print-out  only. 

With  J  down,  a  phase  trajectory  is  printed  for  a  maximum 

of  12  cycles. 
J  -  Set  the  number  of  print-outs  depending  on  zeta.  Change 

the  number  of  overshoot  values  used  to  obtain  the  average 

value. 

JT  -  Use  the  real  time  clock  for  an  automatic  recycle  to  a 
3 

new  parameter  when  have: 

a.  1  l/2  minutes  maximum  after  the  last  print-out. 

b.  30  seconds  maximum  and  ^c  is  1.005  radians  or  less. 


Stop  Switches: 


-S  -  Stops  at  the  end  of  this  parameter  run  or  stops  at  the  end 
of  the  manual  print-out. 
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APPENDIX  B 
00007 

40000 
40001 

40002 
40003 
40004 
40005 
40006 
40007 
40010 
40011 
40012 
40013 
40014 
40015 

40016 
40017 
40020 
40021 
40022 
40023 
40024 
40025 
40026 
40027 
40030 
40031 
40032 
40033 
40034 
40035 
40036 
40037 
40040 
40041 
40042 


COMPUTER  PROGRAM 


ORG 

74  0  00070 

75  0  40701 

FAULT     EXF 

SLJ 

ORG 

74  0  00100 

START     EXF 

75  0  4OOO4 
12  1  40623 

LAZYZETA  thi 

20  0  40425 

STA 

50  1  00000 
75  0  40000 

ENI 

SLJ 

75  0  40674 
50  0  00000 

MANSTART  SLJ 

ENI 

12  0  40565 

LDA 

20  0  40563 

STA 

12  0  40425 

NEWZETA   LDA 

22  0  40564 

AJP 

32  0  40665 

FMU 

32  0  40426 

FMU 

20  0  40404 

STA 

12  0  40624 

LDA 

31  0  40425 

FSB 

22  0  40433 

AJP 

22  2  40433 

AJP 

12  0  40402 

LDA 

20  0  40154 

STA 

20  0  40173 
20  0  40207 

STA 

STA 

12  0  40661 

LDA 

20  0  40400 

STA 

50  0  00000 

ENI 

75  2  40437 

ZZ        SLJ 

50  0  00000 

ENI 

12  0  40652 
12  0  40662 

LDA 

STA 
LDA 

20  0  40617 
12  0  40672 

STA 

LDA 

20  0  40620 

STA 

10  0  00010 

ENA 

61  0  40562 

SAL 

75  4  71000 
50  0  00000 

GOAROUND  SLJ 

ENI 

02  0  40427 

02 

04  0  00000 

04 

12  0  40664 

LDA 

20  0  40156 

STA 

10  0  00000 

ENA 

20  0  00000 

STA 

20  0  40162 

STA 

20  0  40161 
20  0  40157 

STA 

STA 

20  0  40155 

STA 

20  0  40377 

STA 

20  0  40615 

STA 

20  0  40351- 

STA 

20  0  40352 

STA 

20  0  40353 

STA 

20  0  40614 

STA 

20  0-40607 

STA 

20  0  40610 

STA 

20  0  40611 

STA 

20  0  40612 

STA 

20  0  40613 

STA 

12  0  40621 

LDA 

20  0  40326 

STA 

20  0  40335 

STA 

12  0  40622 

LDA 

20  0  40544 

STA 

12  0  40253 

LDA 

20  0  40071 

STA 

12  0  40254 
20  0  40072 

-  LDA 

STA 

75  0  40452 

NEWINPUT  SLJ 

50  0  00000 

ENI 

00007 

0  00076 

0  FAULPRIN 
40000 
00100 


NEWZETA-1 
:TAINDX 


Zl 
ZETA 

0  START 

0  MANUAL 

0  0 

0  NEWSTOP 

0  STOPNOW 

0  ZETA 

0  STOPINDX 

0  TWO 

0  OMEGAN 

0  A 

0  ZETAINDX+1 

0  ZETA 

0  OUIKCOMP 

2  QUIKCOMP 

0  POINT01 

0  TABLES+:. 

0  TABLEL+;. 

0  TABLEM+,. 

0  C0UNT10 

0  TENTHSEC 

0  0 

2  BIGZETA 

0  0 

0  STOP20 

0  STOPZETA 

0  C0UNT12 

0  COUNTZET 

0  EIGHT 

0  DIVIDE 

0  10 

0  STOPRINT 

4  DECOF 

0  0 

0  FORZEROS 

0  0 

0  ONE 

0  UDOT 

0  0 

0  0 

0  THETA 

0  THETADOT 

0  U 

0  T 

0  INDEX 

0  INDEXAVE 

0  PRINTBUF+3 

0  PRINTBUF+4 

0  PRINTBUF+5 

0  LIMINDEX 
LIMITBUF+2 
LIMITBUF+3 


LIMITBUF+4 
LIMITBUF+5 
LIMITBUF+6 
PRINTSET 
0  0K2PRINT+1 
0  BOUNPRIN+4 
0  NEWLIMPT 
0  LIMPRINT+1 
0  FLOTDATA 
0  FLOATEST 
0  FLOTDATA+1 
0  FLOATEST+1 
0  INDXVALU 
0  00000 
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40043 
40044 
400^5 
40046 

40047 
40050 

40051 
40052 
40053 
40054 
40055 
40056 
40057 
40060 
40061 
40062 
40063 
40064 
40065 
40066 
40067 
40070 
40071 
40072 
40073 
40074 
40075 
40076 
40077 
40100 
40101 
40102 
40103 
A0104 
40105 
40106 
40107 


01 

?i 
i\ 

06 
72 
75 
01 
06 
10 
20 
20 
75 
00 
00 
75 

?§ 

50 
12 
20 
20 

1? 
20 

33 
20 

11 

20 

20 
33 
20 
20 

\l 

33 
20 

50 

1? 
30 

32 

33 

11 

22 
12 
22 
50 
10 
20 
20 
20 
20 
75 
00 
00 
75 
50 
75 
50 
00 
00 
75 
50 
13 
32 
30 
22 


8  m&8 

0  40513 
0  40071 
4  71000 
0  00000 
0  40404 
0  00001 
0  40046 
a  71000 
0  40412 
0  00000 
0  40510 
4  71000 
0  40420 
0  00000 
0  00000 
0  40160 
0  40163 
4  60200 
0  40153 
0  40164 
4  60201 
0  00000 
0  40310 
0  00000 
0  40155 
0  40174 
0  40210 
0  40156 
0  40175 
0  40411 
0  40211 
0  40157 
0  40176 
0  40200 
0  40411 
0  40212 
0  40214 
0  40162 
0  40201 
0  40411 
0  40215 
0  00000 
0  40162 
0  40664 
0  40424 
0  40161 
0  40404 
40423 
40053 
40161 
40270 
00000 
00000 
40177 
40202 
40213 
40216 
60200 
40172 
40203 
60201 

00000 
60200 
00000 
0  40206 
0  40217 
4  60201 
0  00000 
0  40215 
0  40411 
0  40201 
0  40112 


0 
2 

0 
3 
0 
0 
0 
0 
0 
0 
4 
0 
0 
4 
0 
4 
0 


k3£ 

01 

06 
RAO 

sw 

06 
RAO 

SW 

06 
COMBINED  ENA 
STA 
STA 
SLJ 
0 
0 
"SLJ 
ENI 
SLJ 
ENI 
LDA 
STA 
STA 
IDA 
STA 
FDV 
STA 
LDA 
STA 
STA 
FDV 
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